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ABSTRACT 

We present an analysis of the Mg II AA2796, 2803 and Fe II AA2586, 2600 absorption line profiles 
in coadded spectra of 468 galaxies at 0.7 < z < 1.5. The galaxy sample, drawn from the Team 
Keck Treasury Redshift Survey of the GOODS-N field, has a range in stellar mass (M«) comparable 
to that of the sample at z ~ 1.4 analyzed in a similar manner bv IWeiner et~aT1 (12009 . W09), but 
extends to lower redshifts and has specific star formation rates which are lower by ~ 0.6 dex. We 
identify outflows of cool gas from the Doppler shift of the Mg II absorption lines and find that the 
equivalent width (EW) of absorption due to outflowing gas increases on average with M* and star 
formation rate (SFR). We attribute the large EWs measured in spectra of the more massive, higher- 
SFR galaxies to optically thick absorbing clouds having large velocity widths. The outflows have 
hydrogen column densities -/V(H) > 10 193 cm -2 , and extend to velocities of ~ 500 km s -1 . While 
galaxies with SFR > 10 M Q yr" 1 host strong outflows in both this and the W09 sample, we do not 
detect outflows in lower-SFR (i.e., logM^/M© < 10.5) galaxies at lower redshifts. Using a simple 
galaxy evolution model which assumes exponentially declining SFRs, we infer that strong outflows 
persist in galaxies with logM*/M Q > 10.5 as they age between z — 1.4 and z ~ 1, presumably 
because of their high absolute SFRs. Finally, using high resolution HST/ACS imaging in tandem 
with our spectral analysis, we find evidence for a weak trend (at la significance) of increasing outflow 
absorption strength with increasing galaxy SFR surface density. 

Subject headings: galaxies: absorption lines — galaxies: evolution — galaxies: ISM 



1. INTRODUCTION 

Galactic-scale gaseous outflows are a basic element of 
the process of galaxy evolution, and are observed in 
galaxies in the local Universe and out to z ~ 6 (e.g. , 
Heckman et all 120001 ; iFrve et alJl2002t lAiiki et al.ll2002|: 
Shaplev et all 120031 : Martini 120051 : iRupke et alJ 12005b: 
Veilleux et aLl 120051 : IWeiner et all 120091 ). They are in- 



voked to explain a wide variety of observational and 
theoretical results concerning the evolving stellar and 
gaseous content of dark matter halos. Outflows are a key 
component of the theory of merger-driven galaxy evolu- 
tion, in which the primary mechanism driving the ob- 
served increase in the number density of "red and dead" 
galaxies from z ~ 1 to today is the merging of gas- 
rich blue galaxies (e.g..lFaber et al. 2007; Tremonti et al.l 
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of merging, and is a favored means by which star for- 
mation is subsequently quenched in the remnant. Out- 
flows must be incorporated into models of disk galaxy 
formation which correctly reproduce the observed scal- 
ing relations between d isk rotation velocity, stellar mas s 
(M*), and radius (e.g.. iDutton fe van den BoscH 120091) . 
Finally, galactic winds may give rise to the Mg II- 
absorbing gas in the extended gaseous halos around 
galaxies observed along QSO sightli nes. Severa l lines 
of evidence suppor t this idea (e.g.. iBond et alJ 120011 : 
IBouche et all I2006L [20071: IMenard fe Chelouchd 12009): 
however, alternative origins for this gas have also been 
suggested, i.e., multiphase cooling of hot halo gas, accre- 
tion of col d gas from the intergalactic medium, or tidal 
stripping ([Mailer fc Bullocld 12001 iTinker fc Ch^ 12001 
IWa^ll99"3j )~ 

In spite of their importance, the physical processes re- 
sponsible for driving galactic outflows are not well under- 
stood. Energy from supernova explosions or AGN feed- 
back heats the surrounding gas and may displace both 
hot and cool gas originating in the interstellar medium 
(ISM), possibly removing it to the galactic halo or be- 
yond. Momentum deposition from radiation or cosmic 
ray pressure may also contribute to the acceleration of 
cool ISM gas; the relat ive importance of t hese two pro- 
cesses is hotly debated (iDekel fe Silk 198 61: iMurray et all 
12003 : ISocrates et aHl2008tlStrickland fc Heckmaii 2009) . 
In galaxies which are known to host outflows in the lo- 
cal Universe, the hot phase is observed in X-ray emis- 
sion, while the cooler phase is detected via optical emis- 
sion lines (e.g., Ha, [N II] AA6548,6583) and in ab- 
sorption against the stellar continuum. However, it re- 
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mains difficult to predict whether galactic winds will 
form and whether they will affect the kinematics of the 
ISM in a given galaxy. Theoretical studies propose that 
there is a "threshold" star formation rate surface density 
(Ssfr) below which supernova-driven superbubbles can- 
not break out of a galac tic disk and form a wind (e.g., 
iMcKee fe Qstrikerl I1977T) ■ Observational constraints on 
the value of a universal threshold Ssfr remain merely 
suggestive (e.g. . iLehnert fe Heckmaij|199d lMartinlll999l ; 
lHeckmanll200a iDahlem et al.ll2006l ) 

For instance, several studies have identified outflows in 
galaxies in which absorption line transitions tracing cool 
gas are blueshifted with respect to the systemic velocity. 
The Na I D AA 5890, 5896 doublet traces the kinematics 
and column of gas at T ~ 10 — 1000 K, revealing out- 
flows in local dwarf starbursts (jSchwartz fe Martinll2004( ) 
and luminous infrared galaxies (LIRGs) out to z ~ 0. 5 
(jHeckman et al.|[2000L lMartirjl2005t iRupke et al.ll2005b[ ) . 
UV spectroscopy of both low and high-ionization transi- 
tions such as Si II A1260 and C IV AA1548, 1550 in Lyman 
Break Galaxies (LBGs) at z ~ 3 has uncovered outflows 
of hundreds of km s _1 (e.g., Sha plev et al.ll2003l ) in these 
objects. All of these galaxies have high spatial concentra- 
tions of star form ation, with Ssfr > 0.1 M© yr _1 kpc -2 
(jHeckmanl 120021 ) : the nearby galaxies have Ssfr values 
among the largest in the local Universe. If there is a 
"threshold" Ssfr required to drive outflows, this sug- 
gests that it is likely equal to or below 0.1 M Q yr _1 kpc -2 
for local galaxies. The value of a threshold Ssfr for 
galaxies in the distant Universe remains poorly con- 
strained, althou g;h the high Ssfr (> 1 Mm yr" 1 kpc" 2 ) 
values in LBGs (jSteidel et alJfl99a iMeurer et al.lll997t i 
suggest a similarly low threshold at z ~ 3. 

Not only are outflows expected to contribute to in- 
creased numbers of post-merger, "quenched" red galax- 
ies, but they may also influence the gradual decline 
of star formation s i nce z ~ 1 am ong blue galaxies 
(jNoeske et al.ll2007at iSato et al.ll200l . as this decline is 
likely driven by gas exhaustion. The presence of winds 
in both distant star-form ing and red sequence galaxies 
has b een established fe.g. JSato et al. 2009; Werner et alj 
[2009J, hereafter W09) , although the co-evolution of out- 
flows with the buildup of M* and the de creasing cosmic 
star formation rate (SFR) density (e.g., Ho pkins! 12004( 1 
remains unexplored at z > 0.5. 

We use a sample of 468 galaxies at 0.7 < z < 1.5 
drawn f rom the Team Ke ck Treasury Redshift Survey 
fTKRS: lwTr~th et alJl200l of the GOODS-N field to ex- 
amine the kinematics of cool (T < 10 4 K) gas traced 
by Mg II and Fe II absorption in coadded spectra. Be- 
cause LIRG-like levels of star formation occur in most 
massiv e star-forming galaxies at z ~ 1 (Noesk e et al.l 
l2007a|) . and because cool outflows are common among 
LIRGs in the local Universe, the galaxies i n our sample 
are pr i me candidates for ho sting outflows (|Rupke et al.l 
120051* iHeckman et al.ll2000f ). Indeed, W09 showed via 
analysis of coadded spectra of ~ 1400 blue, star-forming 
galaxies at z ~ 1.4 that outflows traced by Mg II ab- 
sorption are very common among the objects in their 
sample, and further demonstrated that outflow veloc- 
ities and absorption strengths increase with both M* 
and SFR over the ranges 9.5 < logM^/M© < 10.8 and 
10 M Q yr^ 1 < SFR < 50 M Q yr" 1 . 



However, additional evidence indicates that the out- 
flows observed in galaxies at z ~ 1.4 by W09 must cease 
by z ~ 0.5. W09 suggest that many of the galaxies in 
their sample will ev olve into massive spirals (rather tha n 
ellipticals) at z ~ (|BlantorJl20M iNoeske et al.ll2007bf ). 
The study of Na I kinematics in both star-forming and 
quiescent, r ed-seq uence galaxies at 0.11 < z < 0.54 by 
ISato et all (|2009| ) showed that outflows with velocities 
> 50 km s _1 traced by Na I are detected in ~ 80% 
of the blue cloud galaxies in their sample, while the re- 
ma ining 20% of thes e galaxies lack outflows. Although 
the ISato et alj (|2009fl sample of blue cloud galaxies is in- 
complete, this suggests that at least some star-forming 
galaxies at z ~ 1.4 cease to drive winds in the redshift 
range 0.5 < z < 1.4. 

We use our sample to explore the dependence of the 
equivalent width (EW) of outflowing gas on galaxy 
M* and SFR and examine the relative importance of 
these two properties in determining outflow absorption 
strength among star-forming objects. We follow the co- 
evolution of outflows, galaxy M* and SFR in the range 
0-7 S z ^ 1-4 and search for the expected decline in out- 
flow EW with increasing galaxy age. Finally, using coad- 
ded spectra in concert with the deep HST/ACS imaging 
available in the GOODS-N field, we study the relation- 
ship between Ssfr and the EW of absorption due to 
outflowing gas for the first time beyond the local Uni- 
verse. 

We describe our sample, spectroscopy and imaging 
data in §2. We show evidence for outflow in a few indi- 
vidual galaxy spectra and in a coadded spectrum in §3, 
and describe the measurements we use to quantify out- 
flow properties in §4. Section 5 describes our derivations 
of galaxy properties, and §6 examines trends in outflow 
properties with galaxy SFR, M», SFR surface density 
(Ssfr) and redshift. Section 7 discusses Fe II absorp- 
tion properties. Sections 8 and 9 contain a discussion of 
these results and our conclusions. We adopt a ACDM 
cosmology with Hq — 70 km s _1 Mpc" 1 , S1m = 0.3, and 
SIa = 0.7. Magnitudes quoted are in the AB system un- 
less otherwise specified, and stellar masses are reported 
in units of M . 

2. DATA, SAMPLE, AND STACKING TECHNIQUE 
2.1. Spectra 

We use publicly available spectra from the TKRS 
(jWirth et alj|2004f ) for analysis of Mg II and Fe II kine- 
matics. In brief, the TKRS is a magnitude-limited spec- 
troscopic survey of galaxies selected to have Rab < 24.4 
in the GOODS-N field. Spectr a were obtained us ing one 
hour exposures with DEIMOS (jFaber et al.ll2003f ) on the 
Keck 2 telescope. The 600 1 mm -1 grating blazed at 
7500 A was used with a 1" wide slit. Spectra have a res- 
olution of (Jinst — 1.4 A (jWeiner et al.l 120061 ). and cover 
~ 4600 - 9800 A with 0.648 A pix" 1 . Approximately 
1440 spectra of galaxies and AGN were obtained in to- 
tal. 

The spectra were reduced using an early version of the 
pipeline develope d by the DEEP2 Redshift Survey Team 
(|Coil et aLll2004l J. Newman et al. 2009, in preparation). 
To determine first pass redshifts, the pipeline calculated 
the best fitting source spectrum at each lag position, 
or redshift, from a linear combination of several tern- 
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plate eigenspectra (Strauss et al. 2002; Glazebrook et 
al. 1998). The template spectra included stars, galaxy 
absorption and emission line spectra, and an AGN spec- 
trum. The 10 best fits (i.e., with the smallest \ 2 values) 
were saved, and members of the TKRS team then visu- 
ally checked the best solution, r eplacing it in some cases, 
and assigned it a quality code. Wirth et alj (|2004D per- 
formed a detailed comparison between TKRS redshift 
determinations and redshifts measured in other surveys, 
and estimated that the uncertainty in TKRS redshifts is 
~ 100 km s" 1 or better. 

We limit our sample to spectra with high redshift qual- 
ity codes, i.e., with confidence level > 90%, and with 
coverage of the Mg II AA2796, 2803 doublet. We then 
visually inspect each of these spectra to check for broad 
Mg II emission lines; we exclude four objects with broad 
emission from our analysis (with object IDs 5609, 3660, 
1488 and 9377). These cuts limit the sample size to 625 
objects. 

As noted in lCowie fc Bargerl (|2008l) . some of the TKRS 
spectra suffer from poor sky subtraction, and spectra of 
faint objects can have negative continua. This implies 
that there are large systematic errors not accounted for 
in the one-dimensional error array for each spectrum. If 
the distribution of the errors in the estimates of the sky 
level were symmetric, the sky subtraction errors would 
not affect our results, as we work primarily with coadded 
spectra. However, negative continua (i.e., oversubtrac- 
tions) occur preferentially when there is excess scattered 
light in a given slit, whereas we do not expect that un- 
dersubtraction results from a specific systematic issue. 
Thus, we expect that overestimation of the sky level is 
more likely to occur than underestimation in general. In- 
cluding the objects with oversubtracted sky in coadded 
spectra will tend to reduce the continuum level of the 
coadd. We exclude the 96 spectra with a negative median 
continuum level measured between rest wavelengths 2790 
A and 2810 A (i.e., 15% of the sample). This method of 
tagging spectra with poor sky subtraction is almost cer- 
tainly not comprehensive; however, it likely removes the 
worst cases from the sample. 

Among the remaining 529 spectra, a few have very poor 
wavelength solutions in the blue due to the paucity of arc 
lamp lines at A < 5000 A. This may affect the accuracy 
of the redshift determination as well as the offset of UV 
absorption lines from systemic velocity. In order to assess 
the quality of the wavelength solutions, we fit a single 
Gaussian to the sky lines at 5197 A and 5200 A, which 
are not resolved in the TKRS spectra. We define v = 
km s _1 to be at the average of the true wavelengths 
of the two lines (^5197,5200)) and measure the velocity 
offset of the sky line complex in each spectrum. Our 
assumption that the centroid of a single Gaussian fit to 
the blended lines will be at the average of their central 
wavelengths holds if the lines have equal strengths. In 
the extreme case in which one of the lines is completely 
absent, an additional offset of ±68 km s _1 must be added 
to the calculated offset. To find the mean velocity offset 
of the sky line complex for the sample ((u 6s, 5197+5200}), 
we first exclude spectra with sky lines offset by more 
than 200 km s _1 from ^5197,5200)1 and then calculate the 
mean of the velocity offset in the remaining spectra. The 
resulting (v f> s , 5197+5200} is 13 km s . We then exclude 



all spectra from the sample (61 objects) with sky lines 
offset by more than 130 km s _1 from (v & Sj 5i97+52oo}- 
Accounting for variations in the relative strengths of the 
sky lines, this means that all remaining 468 spectra have 
sky lines offset by no more than 211 km s _1 from their 
true central wavelength; these spectra make up our final 
sample. Figure Q] shows the redshift distribution of the 
objects in the sample. The median redshift of the total 
sample is z = 0.94. 

The standard deviation of the distribution in sky line 
velocity offsets for spectra remaining in the sample is 
40 km s _1 (not accounting for variations in the strengths 
of the sky lines). This dispersion measures the typical 
wavelength uncertainty in the TKRS spectra at A ~ 5200 
A. We estimate considerably smaller uncertainties at red- 
der wavelengths, as arc lines are more abundant in this 
range. For instance, a Gaussian fit to the sky line at 5577 
A yields maximum velocity offsets < 30 km s _1 and a 
dispersion of 10 km s^ 1 for the spectra in the sample. 

The errors in the wavelength solution introduces un- 
certainties in the systemic velocity of the UV absorption 
lines in each spectrum. Throughout the paper, we as- 
sume that the systemic velocity of each galaxy is given 
by the TKRS redshift. As redshift determinations are 
most strongly dependent on sections of the spectra red- 
ward of observed wavelength A bs ~ 5577 A in the red- 
shift range of our sample, they should not be affected by 
poor blue wavelength solutions. However, the UV ab- 
sorption lines of interest are blueward of A b s ~ 5577 A 
in the redshift range z < 1 in the case of Mg II and in 
the range z < 1.15 in the case of Fe II. We estimate the 
error in the systemic velocity for these transitions by as- 
suming the wavelength solution is perfect at A Q b s ~ 5577 
A and that it varies linearly to the blue, with absolute 
velocity offsets at A b s ~ 5200 A between 20 km s _1 
and 200 km s _1 . At the very bluest wavelengths, the ve- 
locity offset is < 135 km s^ 1 for sky line velocity offsets 
< 40 km s _1 , and have a maximum of 675 km s -1 for sky 
line velocity offsets of 200 km s _1 . Near the median red- 
shift of the sample (z = 0.94), the velocity offset at Mg II 
is ~ 27 km s _1 for sky line velocity offsets of 40 km s _1 . 
These errors in the systemic velocities of UV absorption 
lines must be considered when interpreting results from 
coadded spectra, as discussed in §2.31 

2.2. Imaging 

We use the high qu ality HST/ACS imagin g available in 
the GOODS-N field (|Giavalisco et al.ll2004D . The imag- 
ing covers a 10' x 16' area with the ACS F435W, F606W, 
F814W and F850LP bands (B 435 , V 606 , «775 and z 850 ). 
The limiting surface brightness at la in a 1 sq." aper- 
ture in the F850LP ba nd is fiAB = 27.3 mag arcsec -2 
(jGiavalisco et al.1 12004L version 1 release). We use the 
mosaic data in each band with a pixel scale 0.03" pix -1 . 

2.3. Stacking Technique 

We use code written by J.A.N, to coadd TKRS spec- 
tra. The code first masks out bad pixels in each object 
spectrum. It then renormalizes each inverse variance ar- 
ray so that it has a median equal to 1. Each spectrum 
and its associated inverse variance array is linearly in- 
terpolated onto a grid of rest-frame wavelength, and the 
spectra are coadded. The flux in each pixel is weighted 
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Fig. 1. — Histogram showing redshift distribution of all 468 galax- 
ies included in our analysis in bins of Az = 0.05. Only spectra with 
coverage of Mg II AA2796, 2803 are included, which places a lower 
bound on the redshift distribution at 2 fa 0.7. 



by the renormalized inverse variance, so that pixels with 
more noise from sky emission are given less weight, while 
each spectrum overall contributes to the stack in propor- 
tion to its flux and it does not have an extra weighting 
corresponding to its S/N. This results in a coadd which 
is "light- weighted" . This method of coaddition is the 
same as that used in W09. 

Errors in TKRS redshifts as well as errors in the wave- 
length solution will have the effect of broadening the ab- 
sorption and emission features located at the true sys- 
temic velocity in the coadded spectra. The redshift un- 
certainties discussed above are < 100 km s , while ve- 
locity offsets due to poor wavelength solutions are esti- 
mated to be < 135 km s _1 for the 68% of the spectra 
with the smallest sky line velocity offsets, making the 
extreme assumption that the UV absorption lines of in- 
terest are always at the bluest end of the wavelength 
coverage. We therefore expect that absorption observed 
at velocities offset from systemic by > 200 km s _1 in 
both the Mg II and Fe II transitions in the coadd is most 
likely due to absorption at velocities offset from the true 
systemic velocity of the galaxies in the sample. 

3. REST-FRAME UV METAL-LINE ABSORPTION 
IN TKRS GALAXIES 

3.1. Detection of Outflowing Gas in TKRS2158 
The continuum S/N of a typical TKRS spectrum near 
Mg II is ~ 1 pixel -1 . We are therefore not able to vi- 
sually identify Mg II absorption in the vast majority of 
the individual spectra. There are a few brighter galax- 
ies, however, for which rest-frame UV absorption lines 
are evident. 

One example is shown in Figure which includes 
a color HST/ACS image and the full spectrum of the 
object TKRS2158. This object has a redshift z = 
1.43339, has absolute £?-band magnitude Mb = —22.99, 



and has an infrared luminosity logLm/LQ = 11.9 
(|Melbourne et al.l |2005) . Figure [3] includes sections of 
the same spectrum showing UV absorption line profiles. 
The systemic velocity is given by the redshift of the 
galaxy reported in the TKRS, which was determined 
from the velocity of the [O II] doublet. Not only are 
Mg II AA2796, 2803 and Fe II AA2586, 2600 absorption 
lines detected, but Fe II AA2344, 2374, 2382 absorption is 
evident further to the blue. We calculate velocity offsets 
of each line with respect to systemic by fitting a Gaus- 
sian to each and finding the offset of its centroid. These 
measurements are reported in Table[TJ All lines are offset 
from the systemic velocity by ~ —250 km s , with the 
exception of the Mg II 2796 line and the Fe II 2382 line, 
which are offset — 314±19 km s -1 and — 415±36 km s -1 
respectively. These kinematics are a clear indication of 
cool outflow from this galaxy. Discrepancies in the mea- 
surements of the velocity offsets could be due to noise in 
the spectrum such as poorly subtracted sky emission at 
the locations of certain absorption lines. Additionally, 
narrow emission in the Mg II 2796 transition or in Fe II* 
fine-structure transitions may be filling in the some of 
the absor ption profiles and shifting their centers to the 
blue fe.g.. lRubin eTafl 120091 ). 

Figured] shows the spectral regions around the UV ab- 
sorption lines for several more individual TKRS objects. 
Blueshifted absorption lines are evident in a few of these 
spectra, hinting at the presence of outflows. However, 
most of the spectra are not of sufficiently high quality 
to reliably measure velocity offsets or line strengths for 
individual galaxies. 

3.2. Outflow Signature in Coadded Spectra 

Figure [5] shows sections of the coadd of the full galaxy 
sample (468 objects) surrounding the Mg II and Mg I 
A2852 absorption lines. The spectrum has been normal- 
ized by a linear fit to the continuum flux in the con- 
tinua regions —2400 km s -1 < v < —1600 km s -1 and 
800 km s -1 < v < 1600 km s _1 , where v = km s" 1 at 
2803.531 A and 2852.964 A for Mg II and Mg I respec- 
tively. This continuum normalization is applied to all 
coadds presented. The Mg II line profiles in this figure 
are asymmetric, with more absorption on the blue side of 
v = km s _1 . The Mg I absorption line is asymmetric 
in the same sense. Quantitative analysis of these lines 
will be performed in later sections. 

To test whether the properties of the coadd reflect the 
properties of most of the individual galaxies in the sam- 
ple, rather than a few of the brightest galaxies, we mea- 
sure the flux in various velocity ranges with respect to 
Mg II A2803 in individual spectra. This was done in 
W09, and our velocity windows match those used in that 
work: 400 km s _1 < v < 800 km s _1 (window 1), and 
— 150 km s" 1 < v < —50 km s _1 (window 2). Window 

1 provides a measurement of the continuum and window 

2 probes the flux in the deepest part of the blueshifted 
2803 A line (although this velocity range may not probe 
blueshifted absorption in all spectra, due to errors in sys- 
temic velocity determinations as discussed in §2.3|) . We 
plot the average number of counts per pixel weighted by 
the inverse variance in window 2 vs. window 1 in the 
top panel of Figure [6] The points are offset below a 
one-to-one ratio, showing that many of the spectra have 
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Fig. 2. — HST/ACS color image of TKRS2158 in the V606> *775 an d ^850 bands (left) and the TKRS spectrum inverse-variance smoothed 
using a window of 9 pixels (right). The error in the smoothed spectrum is shown in green. The absorption feature at ~ 7650 A is due to 
the atmosphere. Absorption due to Mg II and Fe II is evident at observed wavelengths ~ 5700 — 6900 A. The mean S/N near the Mg II 
transition is 3.4 pix — 1 , The object has an irregular morphology with several bright blue knots, suggestive of a recent or ongoing merger. 
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Fig. 3. — UV absorption lines in the spectrum of TKRS2158. 
The spectrum has been inverse-variance smoothed using a window 
of 9 pixels. The systemic velocity is marked with vertical dotted 
lines, and the dashed lines mark the continuum level. All absorp- 
tion lines shown are blueshifted by at least —250 km s — 1 . Narrow 
redshifted emission is evident in the Mgll transition, resulting in a 
P Cygni-like line profile. 

a count decrement in window 2 (where blueshifted ab- 
sorption is expected) relative to window 1 (where we are 
making a measurement of the counts in the continuum) . 
Approximately 2/3 of the spectra have average counts 
lower in window 2 than in window 1, which indicates 
that a majority of the spectra contribute to the coadded 
flux decrement. 

We also wish to characterize the presence of emis- 
sion in the Mg II transition in individual galaxy spec- 
tra. This is motivated by the work of W09, who iden- 



tify Mg II emission in coadded spectra as well as in 
~ 4% of the individual galaxy spectra in their sample 
at z ~ 1.4. This emission may be due to AGN activity, 
but could also be related to some other physical pro- 
cess. W09 attempt to exclude galaxies which exhibit 
Mg II emission from their analysis of outflow properties 
in coadded spectra, as the complicated continuum shape 
of these individual spectra makes characterization of the 
Mg II absorption line profiles in the coadds difficult. To 
identify individual spectra with Mg II emission in our 
sample, we measure the flux level in the velocity range 
—700 km s _1 < v < —600 km s _1 (window 3) and com- 
pare it to our measurements in window 1 in the bottom 
panel of Figure [5J The placement of window 3 is imme- 
diately to the red of the 2796 A line, and samples the 
region between the absorption lines where Mg II emis- 
sion is likely strongest. Points toward the upper part of 
this plot have "excess emission" in window 3. There are 
7 spectra in our sample with Counts/Pixel in window 3 
greater than 120. Four of these spectra have P Cygni- 
like line profiles at Mg II; the three remaining spectra 
are pushed into the "excess emission" regime because of 
noise in the window. Removal of these spectra from our 
sample does not significantly affect the Mg II line profiles 
in the coadded spectrum. Their inclusion will therefore 
not demand careful modeling of emission in the contin- 
uum, and we retain them in the following analysis. 

4. ANALYSIS OF AVERAGE OUTFLOW 
PROPERTIES TRACED BY MG II ABSORPTION 

In this section we introduce two measurements which 
will be used in the remainder of the paper to assess the 
strength of the absorption due to outflowing gas in coad- 
ded spectra. With both methods we attempt to esti- 
mate the strength of absorption at the systemic velocity 
through analysis of the red side of the Mg II 2803 line. 
We then "correct" the measurement of the absorption 
strength on the blue side of the Mg II 2796 line accord- 
ingly. This corrected measurement should then depend 
on the absorption strength of outflowing gas only. The 
first method uses boxcar EWs for the measurements in 
each velocity range and is useful in the case of low S/N 
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Fig. 4. — UV absorption lines in TKRS spectra, inverse- variance smoothed using a 9 p ixel window and n ormalized to the continuum level. 
Objects are arranged from top to bottom in order of increasing apparent R magnitude (Wirth et al. 2004), and the TKRS ID numbers are 
given in the right-hand panel. The locations of the Fe II AA2586, 2600 and Mg II AA2796, 2803 transitions at the systemic velocity of each 
galaxy are marked with dashed lines. The line profiles in the spectra of objects 8400, 10724, 3230, and 1890 are blueshifted with respect 
to systemic velocity, suggesting the presence of outflows. Dotted lines show the error in each pixel. 
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Fig. 5. — Section of the coadd of all TKRS spectra around Mg II 
and Mg I. The coadd has been normalized to the level of the con- 
tinuum surrounding the absorption lines as described in the text 
and is plotted in black. The dotted line is the error array for the 
coadd, and the black vertical lines mark the systemic velocity for 
each transition. The red line overplotted on the coadd shows the 
Gaussian fits to the red sides of the 2803 A and 2852 A lines; the 
Gaussian model profile for the systemic absorption in the 2803 A 
transition is assumed to describe the systemic absorption in the 
2796 A transition as well. The green spectrum shows the absorp- 
tion remaining after the models have been divided out of the coadd. 
The black horizontal line marks the continuum level. 

coadds, while the second includes more detailed fitting 
of the line profile. The latter method, however, requires 
that the red side of the Mg II 2803 line can be well- 
characterized by a Gaussian. 

Analysis of the absorption strength of outflowing gas is 
complicated by absorption from stellar atmospheres and 
the ISM near the systemic velocity. Photospheric Mg II 
absorption at the systemic velocit y is strongest in F8 - 
Gl type stars (jKinnev et al.lll993l) . in which the EW of 
the 2796 A line can reach ~ 7 A (see |J7]). As outlined 
in W09, very bright late B, A and F stars can exhibit 
asymm etric or shifted M g II absorption due to stellar 
winds (|Snow et al.l 11994 ); these stellar absorption lines 
may be blueshifted by ~ —100 km s -1 . There may also 
be absorption due to the ISM of the galaxies; this gas is 
in the disks of the galaxies and will rarely be blueshifted 
by more than a few tens of km s" 1 (i.e., in the case of a 
rotating disk). The situation for Mg I is simpler; stellar 
Mg I absorption is not blueshifted and is only found in 
the photosphere of F stars. It may be present in the ISM, 
however, and in outflowing gas. 

4.1. Boxcar Method 

In order to characterize the amount of absorption due 
to outflowing gas in each coadd, we make the assump- 
tion that the Mg II doublet is saturated in our spectra, 
such that the 2796 A and 2803 A lines have the same 
depths. Mg II becomes saturated at column densities 
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Fig. 6. — Top: Average counts/pixel in window 2 vs. window 1. 
The points are offset to below a one-to-one ratio (shown in red), 
which indicates that a flux decrement to the blue of the 2803 A 
line is a property of many of the galaxy spectra in the sample. 
Bottom: Average counts/pixel in window 3 vs. window 1. Only 
a few objects are separated from the main locus with high values 
in window 3, indicating there are only a few objects with "excess 
emission". The location of the spectrum of TKRS2158 is marked 
in cyan in both panels. 

(N) of « 10 14 cm -2 , which occur at relatively low hydro- 
gen column densities of > 10 19 cm -2 at solar abundance. 
The ISM and stellar atmospheres typically have columns 
exceeding this value. See §8.1l for further discussion. We 
measure Wdiff, where 
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gas that is not outflowing, and we subtract it from 
W ° to avoid overestimating the outflow absorp- 
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tion strength from the inclusion of absorption due to gas 
at v w km s _1 associated with the ISM and/or stellar 
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atmospheres. 

4.2. Decomposition Method 

In addition to the boxcar measurement discussed 
above, we also adopt the method of W09 to first remove 
the stellar and "stationary" ISM absorption from the line 
profile before making measurements of any outflow. We 
use the model presented in W09: 

F obs {X)=C{l-A sym ){l-A flow ) (4) 

A 8 y m = A279qG(v, X2796, o) 

+ ^2803G(v,A 2 803,O') (5) 

where F a b s (A) is the observed flux density, C is the galaxy 
continuum emission, and A sym and Afi ow are the line 
profiles of the symmetric and blucshiftcd (outflow) ab- 
sorption. Here we make no attempt to include emis- 
sion in Mg II above the continuum in our model; this 
issue will be addressed in more detail in §4.41 A sym is 
made up of the sum of two Gaussians centered at the 
rest velocity of each line in the doublet. In order to cal- 
culate A sym , we fit a Gaussian profile to the red side 
(0 km s- 1 < v < 1600 kms" 1 ) of the 2803 A and 2852 A 
lines. See Figure[5]for a demonstration of this procedure. 
Because the Mg II doublet is blended in our spectra, we 
simply impose the Gaussian fitted to the 2803 A line onto 
the profile of the 2796 A line. The depths of these two 
Gaussians are kept the same (^2796 = ^4-2803 )> since the 
Mg II lines are mostly saturated. It may be that the 
2796 A line is in reality slightly deeper than the 2803 A 
line; in this case we will underestimate the strength of 
the systemic absorption in this line. We divide the coadd 
by this model, and the resulting blueshifted (or outflow) 
absorption line profile is plotted in Figure [5] in green. 

We presume that this profile contains absorption only 
from gas that is outflowing. However, if the doublet ratio 
for the stationary absorption is larger than one, we will 
measure too much outflowing gas (and the emission to 
the red of the 2796 A line will be artificially weakened). 
Redshifted emission in the 2803 A line also affects the 
amplitude and width of the Gaussian we fit to the ab- 
sorption line profile, and can cause an overestimate of 
the EW of outflowing gas. Emission in the 2796 A line 
(evident in Figure [5]) may in turn cause us to underes- 
timate the EW of the outflow (see §4.41 for more discus- 
sion of these effects). Finally, we are assuming that little 
of the absorbing gas is flowing into the galaxies; if we 
were to detect a substantial inflow, we would overesti- 
mate the amount of absorption at the systemic velocity 
due to deeper profiles on the red sides of the lines, and 
underestimate the EW of outflowing gas. 

As can be seen in Figure [5j no absorption is evident 
at positive relative velocities in the line profile shown 
in green, i.e., the "outflow" profile. Some emission oc- 
curs on the red side of the 2796 A line; this will be 
discussed in N4.4I The EW of the feature on the blue 
side of the 2796 A line is 0.51 ± 0.14 A (measured be- 
tween — 1132 km s" 1 < v < km s _1 using a boxcar 
sum). Although the green profile for Mg I is suggestive 
of an outflow, the outflow absorption is not detected, 
with a 3cr upper limit of ~ 0.3 A, possibly indicating 
that the outflow is dominated by gas densities lower than 



those needed to co ntain a significant column in Mg I 
(|Murrav et al.ll2007f) . 

To compute the error on the outflow EW, we first gen- 
erate 1000 realizations of the symmetric absorption pro- 
file by adding noise to the best-fit symmetric model gen- 
erated in the above procedure. The level of the added 
noise is determined by the noise in the coadded spec- 
trum. Each of these realizations is fit with a Gaussian, 
and a double Gaussian profile (A sym ) is created as above. 
We then calculate the standard deviation of the values 
of these 1000 different models at each pixel, which is an 
estimate of the error introduced by the model fitting pro- 
cedure as a function of velocity. This error is combined 
in quadrature with the error in the coadded spectrum 
itself, producing an error at each pixel of the outflow 
profile (Afi ow ). It is then straightforward to use this er- 
ror to calculate the uncertainty in EW measurements of 
the outflow profile. 

4.3. Sensitivity ofWdiff and Outflow EW to Winds 

We now explore the extent to which these two quan- 
tification methods are sensitive to outflows with a range 
of physical parameters. To do this, we generate a series 
of model Mg II absorption line profiles, each with vary- 
ing amounts of absorption at systemic velocity and offset 
to negative velocities. We model both the systemic ab- 
sorption and the outflowing absorption as single velocity 
components with Gaussian optical depth (r) . While cool 
outflowing gas likely consi sts of multiple absorbing clouds 
at different velocities (e.g.. lMartinl2005T ). we are coadding 
our data, and thus we expect that such features will be 
complete ly smoothed out. Com ponents are parametrized 
following fRupke et all (|2005bf ). with variable Mg II col- 
umn density (iV(Mg II)), covering fraction (Cf), Doppler 
parameter (&d), and central wavelength (Ao). In all of 
our models, we choose iV(Mg II) = 10 14 9 cm -2 for the 
systemic component, such that the profile is completely 
saturated, and iV(Mg II)= 10 14 cm~ 2 for the outflow- 
ing component. All components (systemic and outflow- 
ing) have Cf — 0.5, close to the 55% absorption depth 
of the saturated Mg II profiles in W09 (which in coad- 
ded spectra corresponds to the detection frequency of 
outflows times the Cf of cool outflowing clouds). Out- 
flowing components are given relative velocities of -100, 
-200, and -300 km s _1 ; we also create models with no 
outflowing component. In each model, the systemic and 
outflowing components have the same bo] however, we 
allow this parameter to have the values bp = 50, 150, 
and 250 km s _1 in different models. We therefore have 
a grid of 3 x 4 models in bo and outflow velocity space. 

We smooth each of these models to the velocity 
resolution of the i ndividual spectra, Cmst = 1.4 A 
(jWeiner et al.l l2006h . adjusted to the rest-frame at the 
median redshift of the sample (z = 0.94). A velocity 
resolution element in the coadded spectra is larger than 
this, due to uncertainties in the redshift determinations 
and wavelength solutions for the individual spectra (see 
§2.ip ; however, the results are similar if we repeat the 
analysis using a velocity resolution that is 2a ms t- After 
smoothing we add different levels of noise such that the 
resulting spectra have S/N = 3, 6, and 9 pix -1 , con- 
sistent with the range in S/N levels of coadded spectra 
we create in Sj6] (see Table [3]). We generate 1000 realiza- 
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tions of each model at each S/N level, and then measure 
Wdiff and perform our decomposition analysis for each 
realization. 

Figure [7] shows the distribution in the measured out- 
flow EW and Wdi / / values for all models as a function of 
the model outflow velocity. The points have been offset 
by an arbitrary amount in velocity so that they do not 
overlap. The mean outflow EW and Wdiff for models 
with S/N = 3 pix -1 , 6 pix -1 and 9 pix -1 are shown 
with cyan diamonds, blue circles, and black squares re- 
spectively. The point size increases with bo- The error 
bars show the 90% confidence intervals in the measured 
quantities. The size of these intervals and the central val- 
ues indicate the degree to which outflow EW and Wdi f f 
are successful in characterizing the underlying physical 
absorption profile. For instance, at S/N = 3 pix -1 , 
Wdi f f and outflow EW have very broad distributions for 
all models. For models with no outflow, both of these 
measures exceed values of 0.5 A in between 12% and 
40% of realizations, with outflow EWs for models with 
bjj = 150 km s -1 and 250 km s -1 exceeding ~ 1 A in 
6% and 12% of realizations, respectively. Likewise, mod- 
els with outflows can have Wdiff or outflow EW < 0.5 
A, although in general the distributions shift to higher 
values of these quantities as outflow velocities increase. 
We conclude, however, that in the case of coadds with 
S/N ~ 3 pix" 1 , we are able to confirm the presence of 
outflows only if the measured values of outflow EW are 
> 1.2 — 1.3 A and Wdiff values are > 0.8 A. These values 
are recovered for outflow velocities > 100 km s -1 , and 
occur more frequently with higher bo- If lower values 
are measured, the presence or lack of outflows remains 
ambiguous. 

These measurements become more sensitive with in- 
creasing S/N. At S/N = 6 pix -1 , Wdiff and outflow 
EW can be as high as 0.5 - 0.9 A for models with no out- 
flow. However, models with outflows yield mean values 
higher than these if bp = 150 or 250 km s -1 . Values of 
outflow EW > 0.85 A are measured for at least 95% of 
realizations of models with high values of bo and out- 
flow velocity. At S/N = 9 pix" 1 , models with no outflow 
have Wdiff and outflow EW < 0.4 — 0.7 A, while models 
with outflows and bu = 150 or 250 km s -1 again yield 
higher values in at least 92% of realizations. Models with 
6_d = 50 km s -1 have higher mean values if the outflow 
velocity > 200 km s -1 . 

Although this modeling relies on a very simplistic pa- 
rameterization of the Mg II absorption profiles in our 
coadded spectra, it suggests that we are sensitive to sat- 
urated Mg II-absorbing outflows with velocities on the 
order of 100 km s -1 in coadds with S/N ~ 6 — 9 pix -1 . 
We are sensitive to outflows in S/N ~ 3 pix -1 coadds 
only with large velocities and bjj . These findings will be 
discussed further in [JS 

4.4. Mg II in Emission 

We observe emission in the Mg II transition in both 
individual galaxies (see the P-Cygni profile in Figure [3]) 
and in coadded spectra (e.g., Figured]). Emission in the 
latter is obvious after the decomposition of systemic and 
outflow profiles is performed and one can identify the 
decrement of absorption on the red side of the 2796 A 
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Fig. 7. — Top: Distribution in measured outflow EW values for 
model spectra as a function of the input outflow velocity. The 
points are offset slightly in velocity to prevent overlap. The mean 
outflow EW for realizations with S/N = 3 pix -1 , 6 pix -1 and 
9 pix -1 are shown with cyan diamonds, blue circles, and black 
squares respectively. Small, medium, and large symbols show re- 
sults for models with bo = 50, 150 and 250 km s — L The error 
bars indicate the 90% confidence intervals for each set of model 
realizations. Bottom: Same as above, for Wdiff- 

line as compared to the red side of the 2803 A transi- 
tion. This emission is observed in z ~ 1.4 star- forming 
galaxies (W09), a l uminous starburst galaxy at z ~ 0.7 
(| Rubin et al.l [2009), as well as in narrow-line Seyfert 
galaxies such as two of th e ultraluminous inf r ared g alax- 
ies (ULIRGs) studied in IMartin fc Bouchl pOOl . A 
similar P-Cygni-like profile in the Na I transition was 
observe d in NGC 180 8, a starburst galaxy driving an 
outflow !|PhilUDslll993T) . While the origin of this emission 
remains unclear, we suggest it may be at least in part 
due to resonance-line scattering off of the receding side 
of an expanding shell related to the o bserved outflow, as 
in the case of Lya emission in LBGs (jPettini et al.ll200ll 
K. Rubin et al., 2009, in preparation). 

To better understand the effect of this emission on 
our outflow EW measurements, we develop a alternative 
method of characterizing the Mg II doublet profiles in 
coadded spectra. We again assume that the absorption 
at the systemic velocity produces a saturated Gaussian 
absorption profile as in Equation [5] We then further 
assume that there is additional emission on top of this 
continuum on the red side of each line with a Gaussian 
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profile. The amplitudes of the emission lines have a ratio 
2:1, and the lines have a variable velocity offset with re- 
spect to systemic. In a true P-Cygni profile, the velocity 
at the peak of the emission (and the shape of the profile 
in general) depends on a number of factors, e.g., the out- 
flow geometry and C / , and the velocity dispersion of the 
outflowing gas. Our model can be written as follows: 

FobsW = C(l — A sym )(l — Afi ow ) + F em (6) 

F em ,W — A 2 79 6 G(v - V , A 2 796, o) 

+ A 2 803G(v - V , A 2 803, o) (7) 

where F ern is emission in excess of the continuum and 
^4-2796 = 2 ^2803- We may then fit this model to the red 
sides of both lines in the doublet simultaneously (where 
A flow = 0), in the velocity ranges km s _1 < v < 
268 km s" 1 and km s" 1 < v < 1600 km s" 1 for the 
2796 A and 2803 A lines, respectively. We subtract the 
fitted emission model (F em ) from the coadd, divide out 
the model symmetric absorption, and measure the boxcar 
EW in the range —1132 km s _1 < v < km s _1 for the 
resulting profile. 

Because of the number of free parameters in this model 
(5), this method does not generally produce acceptable 
fits for the coadded spectra in our study, as the results 
are often driven by noise features in the line profiles. 
However, the model does successfully characterize the 
red sides of the doublet lines in the much higher-S/N 
coadds of W09. We find that in the W09 coadds with the 
strongest emission features (e.g., the lowest- and middlc- 
M* subsamples), the outflow EW we calculate using this 
method is ~ 6 — 17% lower than the outflow EW calcu- 
lated in W09 (and described above in §4.2|) . In coadds 
such as the high-M* W09 subsample, the difference in 
outflow EWs calculated with the two methods is < 3%. 
While this model is quite simplistic, these results indicate 
that Mg II emission causes an overestimate of the out- 
flow EW calculated using our standard method fi )4.2|) . 
A more complete, physically-motivated model including 
radiative transfer is required to fully quantify this effect, 
and additionally may more tightly constrain other char- 
acteristics (e.g., radial extent, density) of the cool gas 
outflow. 

5. PROPERTIES OF SAMPLE GALAXIES 

We wish to characterize outflow absorption strength 
at 0.7 < z < 1.5 in galaxies with a range in SFR, M* 
and Ssfr as well as explore the evolution of outflows. 
To do this, we create subsamples of the galaxies with 
similar SFRs, M*s, etc., and coadd spectra within a given 
subsample. We then compare outflow properties among 
these coadds. Here we describe how photometry, sizes, 
M», SFR, and quantitative morphologies are derived for 
our sample. 

5.1. Rest-frame Colors and Luminosities 

We use photometry fro m iWeiner et all (I2006T) de- 
rived from ACS i maging (Giavali sco et al.1 I2004D and 
ICapak et al.l (|2004fl ground-based photometry, and con- 
verted to absolute Mb and r est-frame U — B color using 
the K-correction routine of iWillmer et al.l |2006) . Er- 
rors in the observed optical magnitudes and colors are 



0.05-0.07 mag and 0.07-0.1 mag, respectively. The la 
errors introduced by the K-co rrection procedure a re 0.12 
for M B and 0.09 m U — B (jWeiner et al.ll2006fl . The 
left-hand side of Figure [8] shows a color-magnitude di- 
agram (CMP) for o ur sample. The solid line is from 
Will mer et all (|200ll . and marks the division between 
the red sequence, or the narrow region populated by 
early- type E/SOs in the CMD, and the blue cloud, or the 
wider area in the CMD populated by bluer spirals and 
separated from the red sequence by a narrow "valley" in 
the surface density of objects. 

5.2. Sizes 

We wish to measure Ssfr f° r °ur sample, where 
Ssfr = SFR/-7rRg F t , 2 , for comparison with the sug- 
gested local threshold for driving outflows, Ssfr > 
0.1 M Q yr _1 kpc~ 2 . In the local Universe, measurements 
of the SFR from Ha, far-IR or extinction-corrected UV 
fluxes are combined with measurements of the sizes of 
star-forming regio ns (RsFq/2) to determine Ssfr (e.g., 
iMeurer et al.lll997tlLehnert fc Heckmanll 995). The sizes 
are constrained using measurements of the half-light ra- 
dius from Ha or UV (e.g., 2200 A) imaging, which 
are direct tracers of nebular emission or emission from 
young stars. Local star burst galaxies with high val- 
ues of SFR per unit star-forming surface area (several 
M yr _1 kpc -2 ) have values of SFR per unit surface 
area in the optical disk, i.e., SFR/ ttD 2 ,^ 9 , which are at 
least 3 orders of magnitude lower (|Martinl H999). It is 
therefore important to probe the spatial extent of UV 
flux when making size measurements of distant galaxies 
for comparison with local results. Continuum emission 
at redder wavelengths arises from older stellar popula- 
tions, which are not expected to be important for driving 
large-scale outflows. 

We use the half-light radii of t he TKRS galaxies mea- 
sured by IMelbourne et al.l (|2007f) to parametrize galaxy 
size. These authors fitted successively larger elliptical 
apertures to the ACS images for each galaxy and calcu- 
lated the fluxes and intensities within each aperture. An 
iterative curve-of-growth analysis was used to determine 
the flux level of the sky. From these measurements, the 
total flux of the object was calculated. Apparent half- 
light radii are equal to the semimajor axis of the ellipse 
that contained half of the total flux, and were corrected 
for the point-spread function of the image. To determine 
half-light radii at 2200 A in the rest-frame, we interpo- 
late between the radii measured in the passbands to the 
red and blue of 2200 A in the rest-frame of each object, 
assuming that all of the light measured in each band 
is observed at the central wavelength of the filter (4297 
A for the -B435 band and 5907 A for the VI306 band). In 
cases in which 2200 A in the rest-frame is blueward of the 
i?435-band filter, we simply adopt the £?435-band radius. 
We then use the angular diameter distance to compute 
the rest-frame UV half-light radius in kiloparsecs (R1/2). 
IMelbourne et al.l (|2007t ) estimate that they obtain accu- 
rate radii to within < 10%. This level of uncertainty 
applies strictly to the rest-frame i?-band half-light radii 

9 D25 is the apparent isophotal diameter measured at surface 
brightness = 25 mag/D" 
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they derive by combining the radii measured in the ob- 
served bands in a weighted mean, with weights depen- 
dent on the overlap of each observed passband with the 
rest-frame _B-band. We assume this level of uncertainty 
applies to our R4/2 as well. 

We choose to use R1/2 to parametrize the sizes of the 
star-forming regions in our galaxies for its simplicity; 
however, this measure may in fact significantly overesti- 
mate the size scales relevant for driving outflows. Many 
of the galaxies have extended and clumpy morphologies, 
such that R4/2 is quite large (> 10 kpc), while much of 
the UV emission arises in a few small but widely sepa- 
rated bright knots. On the other hand, the distance be- 
tween star-forming knots may be intimately connected 
to the morphology of the outflowin g gas itself (see, e.g., 
iHeckman etHI Il990; Martin 200$ for some discussion 
of outflow morphology). Future studies of outflows in 
individual galaxies will warrant more careful analysis of 
the size scale and distribution of star formation. 

5.3. Stellar Mass 

Near-I R photometry of the GOODS-N field was pub- 
lished bv IBundv et~aLl (|2005f ). who derived M* for 202 
objects in our sample (out of a total of 468 ). As we did 
not wish to limit this study to only objects with K-band 
derived M„, we use a calibration derived in W09 to con- 
vert re st-frame color and magnitudes into M* . iBell et al~l 
(2003) compute the relation between rest-frame color 
and M*/L B using SDSS and 2MASS photometry of local 
galaxies. They give this relation for a "diet Salpeter" 
IMF: 

\ogM*/L B {z = 0) = -0.942 + 1.737 (B - V Vega )- 

This must be adjusted according to the redshift of each 
galaxy. W09 derive a redshift correction to this relation 
using the if -band magnitudes and M* (with Chabrier 
IMF) a vailable for 11924 objects in the DEE P2 redshift 
survey ([Davis et al.ll2003l : IBundv et al.|[2006h . These ob- 
jects lie in a redshift range < z < 1.5. W09 performed 
a least-squares fit between M* derived from rest-frame 
color and those derived from if-band photometry, and 
give a correction term C'k '■ 

M*=L Byega x M*/L B (z = 0) 
x C K {U-B,z); 
log C K (U -B,z) = -0.0244 - 0.398z 

+ 0.105 {U-B Vega ). 

The authors find a scatter of 0.25 dex about the fit for 
Ck- We apply this correction to our data, and compare 
if -band M* to the corrected color-derived masses where 
possible (for 202 objects). We find there is a 0.09 dex 
mean offset and a dispersion of 0.22 dex between the 
two stellar mass estimates. We use the corrected color- 
derived masses with a Chabrier IMF for the full sample 
in the following analysis. 

5.4. Star Formation Rate 

The rich multi-wavelength data set in the GOODS- 
N field makes several different methods available for the 
derivation of the total SFR for the galaxies in our sample. 
[O II] li ne luminosities hav e been measured for TKRS 
galaxies (jWeiner et al. 2006). 24 /im fluxes and IR-based 



SFR s for a subset of objects in our sample are available 
from lMelbourne et al.l {2005). In order to make compar- 
isons with W09, we adopt the ir method of determ ining 
SFR. The BV photometry of lCapak etUI (|2004h mea- 
sures the flux at 1800 A - 2800 A and 2200 A - 3400 A 
in the rest-frame for our redshift range. We use these 
measurements to derive absolute magnitudes at 1500 A 
and 2200 A u sing t he K- correction code des cribed in 
iWillmer et all (|200ll and iWeiner et all pOOST ). Using 
these luminosities we find the slope of the UV contin- 
uum, 0, and calculate the attenuation fr om this slope 
using the relation Apjjy = 3.16 -I- 1.69/? (jSeibert et al.1 
I2005t iTrever et al.l 120071) . There is a ±0.9 mag scatter 
in this relation, which results in a 0.36 dex uncertainty 
in the SFRs. We use the values of Apuv to calculate 
an un extincted UV lu minosity, and in turn calculate the 
SFR (|Kennicuttlll99l for a Kroupa IMF: 

SFR UV (M yr- 1 ) = 1.0 x 10- 28 L i5ool (erg s^Hz- 1 ). 

To investigate systematic errors in these SFR esti- 
mates, we com pare them to SFRs der ived from the IR 
luminosities of Melbou rne et alj (|2005f ). These authors 
used the publicly available MIPS imaging of the GOODS- 
N field t o measure 24 ^m fluxes to a limit of 25 fiJy, and 
used the iLe Floc'h et alJ (|2005h prescription to convert 
24/zm flux to total IR luminosity (L IR ). 213 of the 407 
blue cloud galaxies i n our sample ha ve Lm estimates 
available. We use the IBell et al.l (|2005f ) relation between 
L IR and SFR with a Kroupa IMF to calculate SFR(IR) 
for these galaxies, ignoring the contribution from unex- 
tincted light from young stars in the UV, which we ex- 
pect to be small (jBell et al.ll2005l W09). We find that 
for IR-selected galaxies, which may have slightly larger 
IR-to-UV emission ratios than is typical (W09), Lj R - 
derived SFRs are 0.21 dex higher than the UV-derived 
SFRs, with a dispersion of 0.43 dex. The expected off- 
set will likely be smaller for galaxies not detected in the 
MIPS imaging. For consistency with W09, and because 
we lack measurements for nearly half of the blue 
cloud galaxies in our sample, we adopt the UV-derived 
SFRs described above in our analysis. However, these 
uncertainties must b e considered when absolute SFRs are 
discussed, as in §6.21 IBell et afl (|2005h report systematic 
and random errors in their IR / UV derived SFRs of 0.3 
and 0.4 dex; we must account for these errors in addition 
to those introduced from using purely UV-derived SFRs. 

The right-hand side of Figure [8] shows a plot of 
logSFRuv vs - logM* for all objects in our sample with 
UV SFR measurements. Upper limits on the SFRs in 
red sequence galaxies are marked in red. While our sam- 
ple includes galaxies in a range of M* comparable to the 
sample of W09, the mean SFR in W09 is ~ 0.57 dex 
higher than in this study. 

5.5. Morphology 

We use the Gini (G) and M20 measurements made for 
TKRS galaxies in the «775-band by J. Lotz (2008, private 
communication) to quantify galaxy morphol ogy. These 
are nonparam etric measurements described in lLotz et alJ 
(|2004 120061) . G quantifies the relative distribution of 
light among a galaxy's pixels, and is high if there are 
only a few very bright pixels. M20 is the second-order 
moment of the brightest 20% of a galaxy's pixels; it is 
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Fig. 8. — Left: Color-magnitude diagram for all objects in our 
sample. The solid line marks the division between the red sequence 
and the bl ue cloud objects (in red and black, respectively); it is 
taken from lWillmer et al.l J20TO . Right: SFR UV vs. M* for the full 
sample. Downward red arrows mark upper limits on the SFRtjv 
for objects on the red sequence. 

larger in galaxies in which the brightest pixels are fur- 
thest from each other. Figure shows the distribution 
of these parameters for our sample (excluding 115 which 
lack high-quality G/M20 measurements). The dividing 
lines for d i fferen t morphological types are taken from 
lLotz et ail (|2008f) . who applied these divisions to galax- 
ies at 0.2 < z < 1.2. There are 63 mergers, 226 late-type 
objects, and 64 early-type objects in our sample. 

6. DIVISION OF SAMPLE BY GALAXY 
PROPERTIES 

6.1. Division by SFRtjv > M* ; and Morphology 

To examine trends in Wdiff (where Waff = 
W —W ) with SFRtjv, M*, and morphol- 

2796A,bluo 2803A,rcd ; uv ' ' 1 

ogy, we have divided our spectra into several different 
subsamples and coadded them. We calculate Wdiff for 
each subsample. These subsamples and measurements 
are listed in Table [H 

First, we coadd only galaxies on the red sequence. The 
coadd has a very low S/N = 1.4 pix -1 in the continuum 
surrounding Mg II, and so cannot be used to examine 
outflows in these objects. Because of their low S/N, we 
exclude all red sequence galaxy spectra from the subsam- 
ples described in the following, with the exception of the 
morphologically divided subsamples. 

We choose to divide the TKRS galaxies by the 25th 
and 75th- percentile values of logSFR (0.555, 0.942) and 
logM, (9.86, 10.49). Coadds of the spectra for these sub- 
samples are shown in Figure [TTJ] Figure [TT] shows Wdi / / 
for the SFRyv-divided and M» -divided subsamples. The 
outflow absorption strength (as quantified by Wdiff) 
rises significantly with SFRtjv and M* between the mid- 
dle and highest SFRtjv an d M» subsamples. There ap- 
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Fig. 9. — G vs. M20 measured in the 1775-band for objects in 
our sample. Measurements were provided by J. Lotz (2008, pri- 
vate communication). The colored lines show the regio ns occupied 
by dif ferent morphological types and are adopted from lTxitz et alj 
(2(508). Red sequence galaxies are marked in red. 



pears to be detected outflowing gas in the lowest-M* sub- 
sample; however, this coadd has S/N ~ 3 pix -1 , and thus 
is unlikely to yield reliable measurements (see N4.3I for a 
discussion of this issue) . We also perform our decomposi- 
tion analysis on each of these coadds, shown in FigurefTUl 
Measurements of outflow EW for each coadd are shown 
in Figure 1111 and show consistency with the Wdi f f mea- 
surements. Outflow EW results from W09 are shown 
in blue. Note that measurements of outflow absorption 
strength are significantly higher for the W09 z ~ 1.4 
galaxies than for the middle-M* TKRS galaxies. As dis- 
cussed in H5.41 while the TKRS and W09 galaxies have 
a similar range in M*, on the whole the TKRS galax- 
ies have a mean SFRtjv that is lower by ~ 0.6 dex. In 
addition, the highest-M* TKRS galaxies include objects 
with the highest SFRs in the sample. This suggests that 
either outflow absorption strength is most closely corre- 
lated with SFR, or that there is evolution of outflows in 
logM* ~ 10 galaxies between z ~ 1.4 and 1. This will 
be discussed in greater detail in fj8] 

The coadd of early-type galaxies (as classified by G- 
M 2 o) has very low S/N and is not useful for measur- 
ing outflow absorption. Wdiff for the merger candidate 
and late-type subsamples are the same within the er- 
rors, at Wdiff ~ 0.3 — 0.5 A; i.e., there does not appear 
to be a significant difference in the strength of outflow 
absorption in late-type galaxies and merger candidates 
(see Tabled]). However, it is difficult to disentangle the 
effects of morphology and SFRtjv ° n outflow strength 
in this analysis. Studies of galaxy mergers (i.e., LIRGs 
and ULIRGs) in the local Un iverse show that they host 
except ionally strong outflows (jMartin 2005; Rupk e et al.l 
2005b) in comparison to local late-type galaxies. On the 
other hand, these merger remnants also have some of the 
highest SFRs at z ~ 0. This degeneracy between mor- 
phology and SFR is broken at z ~ 1, where the majority 
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of LI RGs have disk-like morphologies (Melbourn e et al] 
2008); thus galaxies at z ~ 1 may provide the ideal lab- 
oratory for investigating the effects of these two param- 
eters on outflows. Higher S/N spectra are required to 
examine these effects in greater detail. 

6.2. Division by Ssfr 

To test for a dependence of outflow strength on SFR 
surface density (Ssfr) at z > 0, we assume that 
Rsf,i/2 = R-i/2) an d combine this with our measure 
of SFRuv to calculate the global (i.e., flux averaged) 
Ssfr = SFR/7rRg F j, 2 for each object. Figure [T"2l shows 
the distribution of log SFRuv an d log ttR 2 ^ ^ or ^ ne ^ me 
galaxies in the TKRS sample. The lines show the 25th 
and 75th percentile values of log Ssfr. (-1-36, -0.66), 
which we use to subdivide the sample. Figure [TBI shows 
the coadds of the spectra in these subsamples and the 
results of the decomposition procedure. 

Figure [TTl shows Wdiff and outflow EW vs. Ssfr in 
each coadd. There is a slight increase in outflow absorp- 
tion strength with Ssfr evident between the middle- and 
high-SspR subsamples. The lowest-SsFR subsample has 
S/N ~ 3 pix -1 , and so here we are sensitive only to high 
velocity outflows with large velocity widths. As discussed 
in §4.3[ we can be confident that outflows are present in 
subsamples with low S /N coadds only when the measured 
outflow EW > 1.2 - 1.3 A. The outflow EW we measure 
for the lowest-EsFR bin is only slightly below this limit, 
and is therefore only suggestive of the presence of out- 
flow. Comparing Figures [TTJ and [T4j it appears that the 
outflow absorption strength is higher in the high-SFRuv 
subsample than in the high-SsFR subsample, suggesting 
that outflow absorption strength is more strongly corre- 
lated with absolute SFRuv than with Ssfr- 

We also consider our results in the context of the sug- 
gested local "threshold" for driving outflows, Ssfr ^ 
0.1 M© yr^ 1 kpc~ 2 . Because the existence of and a pre- 
cise value for a strict threshold Ssfr for driving outflows 
have not yet been obs ervationally established (see, e.g., 
iStrickland et al.l r2004f ). it is interesting to search below 
the suggested threshold for evidence of winds. Given 
the uncertainty in our absolute SFR determinations, it 
is difficult to differentiate which of our individual galaxies 
have Ssfr above or below Ssfr = 0.1 M© yr _1 kpc~ 2 . 
The solid green line in Figure [T2] shows a line of con- 
stant Ssfr = 0.1 M Q yr _1 kpc -2 for comparison with 
the distribution of our sample in 7rR 2 y 2 -SFRuv space. 
From the placement of this line, it appears that many 
of the galaxies in the middle-SsFR subsample have Ssfr 
values below the threshold. If we shift this line by -0.21 
dex in SFRuv to refle ct ou r possible systematic under- 
estimate of SFR (see §5.4[) as indicated by the dashed 
green line, a smaller fraction of the sample falls below 
the threshold. A 0.21 dex correction to the Ssfr values 
of the sample is likely too large, as this offset is derived 
from a comparison between IR- and UV-based SFRs for 
the IR-detected galaxies only. However, we apply this 
correction and create a new subsample with "corrected" 
Ssfr < 0.1 M Q yr _1 kpc -2 . The coadd of these spectra 
(referred to as the log Ssfr, Corrected sample in Table [2]) 
yields S/N - 4.7 pix" 1 , W dlff = 0.64 ± 0.24 A and out- 
flow EW = 0.90 ±0.29 A. This implies that some galaxies 



below a flux averaged Ssfr = 0.1 M yr _1 kpc -2 drive 
outflows. Because there is a la dispersion in the com- 
parison between SFRuv and SFR(IR) of 0.43 dex, some 
fraction of galaxies in this subsample do have true Ssfr 
values above the threshold; however, the absorption line 
profile in the coadded spectrum reflects the mean absorp- 
tion properties of the subsample, which has Ssfr below 
the threshold in the mean. 

6.3. Redshift Dependence 

As noted in $5A\ the TKRS and W09 samples span 
a similar range in M», but the TKRS sample is offset 
to lower SFRs by ~ 0.6 dex. To examine evolution in 
outflow properties with redshift, we compare our out- 
flow measurements to those from W09. The 25th and 
75th-percentile M* divisions of our sample are quite close 
to the divisions used in W09 (logM* < 9.88,9.88 < 
logM* < 10.45, and logM* > 10.45). To further inves- 
tigate similarities between the two sets of subsamples, 
we compare the distributions of specific SFR (SFR/M*; 
SSFR). Figure [TJ shows the SSFR distribution of TKRS 
(solid line) and DEEP2 (dotted line) galaxies in each bin. 
Each subsample has a symmetric distribution in SSFR at 
both redshifts; these distributions are simply offset from 
each other because of the decline in global SFR with 
decreasing redshift. The symmetry of the distributions 
suggests that none of the subsamples are severely con- 
taminated with non-star-forming galaxies. See §8.41 and 
Figure [19] for further comparison of the SSFR-M, rela- 
tions of the W09 and TKRS samples and a discussion of 
the evolutionary connection between the galaxies at the 
two redshifts. 

Although the TKRS and W09 samples are nearly dis- 
joint in SFRuv-M* space, we may construct a compari- 
son subsample of 43 TKRS galaxies which lie above the 
lower envelope of the SFRs for the W09 galaxies. Fig- 
ure [TBI shows the TKRS (black circles) sample and a ran- 
dom selection of half of the W09 (green crosses) sample 
in logSFRuv-logM, space. We include all TKRS galax- 
ies above the solid line in a high-SFRuv, "DEEP2-like" 
sample. Figure [17] shows the coadd and decomposition 
analysis of this subsample. 

The red point in Figure [TTJ shows results for this 
DEEP2-like sample. This sample as well as the high- 
SFRuv TKRS subsample have outflow EWs similar to 
or larger than the W09 coadds. This confirms our pre- 
vious finding that galaxies at z ~ 1.4 and 1 with similar 
SFRs have strong outflow absorption, and that galaxies 
with lower SFRs have weaker outflows. We note that 
the galaxies in the highest-M* TKRS subsample have a 
median log SFRuv = 0.99, with 68% of the subsample 
in the range 0.78 < log SFRuv < 1-22. These galaxies 
therefore have SFRuv values in a range similar to the 
W09 sample; the high outflow EW for this subsample 
is consistent with the measurements for our other high- 
SFRuv subsamples. 

7. FE II 2586, 2600 ABSORPTION 

We now focus on measurements of Fe II absorption 
in our spectra. A constraint on the extent to which 
Fe II 2586, 2600 absorption is detected with Mg II in 
outflowing g as is valuable for a num ber of reasons, as 
discussed in iMartin fe Bouchel (|2009fh While Mg II is 
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Fig. 10. — Sections of the coadded spectra for subsamples divided by SFRtjv (left) and M* (right) around Mg II (black). The coadds 
have been normalized to the level of the continuum surrounding the absorption lines as described in the text. The symmetric absorption 
profile is shown in red; the outflow profile is shown in green. The error in each pixel is shown with the dotted lines, and the black vertical 
lines mark the systemic velocity for each transition. The number of spectra in each subsample is shown in the upper right of each panel. 



found in a wide range of gas densities, Fe II is present 
only at higher densities (in more neutral gas) , and so may 
provide information about outflow gas density. Second, 
the Fe II 2586, 2600 transitions have oscillator strengths 
0.11 and 0.39 times that of Mg II 2796, and Fe is less 
abundant in general than Mg; thus Fe II absorption lines 
are generally less strongly saturated than Mg II lines. 
They may therefore be used to place more stringent con- 
straints on the cool outfl ow column density. For instance, 
iMartin fe Bouche! (|2009h find that the Fe II absorption in 
ULIRG outflows at z ~ 0.25 requires optical depths 2-3 
times larger than those derived from analysis of Mg II. 
Finally, because Mg is generated in Type II supernovae 
in younger stellar populations than Fe, comparison of the 
relative abundances of Mg II and Fe II in the outflow may 
constrain the enrichment history of the gas. 

The Fe II 2600 transition is always weaker than Mg II 
2796 in QSO absorption line systems, with at most 2/3 
the EW of Mg II 2796 (|Churchill et alJl2o5ch . This is 
also true in stellar atmospheres. Fe II 2600 reaches its 
maximum absorption strength in F9 - G2 stars, while 
Mg II reaches its m aximum strength in F8 - Gl stars 
(|Kinnev et al.lfl993l although asymmetric Mg II absorp- 
tion is found in B-F stars as mentioned in In early 
type stars, Fe II lines do not exhibit P-Cygni profiles or 
asymmetric absorption; instead if mass-los s effects are 
obser ved the lines are simply blueshifted (|Snow et al.1 
H99l . 

We use UVBLUE theo retical stellar spectra 
([Rodriguez- Merino et al.l l2005h to measure the relative 



EWs of these transitions in stellar atmospheres. We first 
normalize these spectra, calculating the continuum in the 
same regions used for our galaxy coadds. We measure 
EWs between -500 km s" 1 < v < 500 km s _1 for Fe II 
2600 and between -500 km s _1 < v < 385 km s" 1 for 
Mg II 2796. Note that these spectral regions contain ab- 
sorption not only from the named transitions, but from 
several weaker transitions also found in stellar atmo- 
spheres (the Fe II* 2599 transition, for example). Using 
solar metallicity models, we obtain EW(2796) ~ 6.8 A 
and EW(2600) ~ 3.1 A for a solar type star, yielding 
a ratio of ~ 2.2. As the atmospheres become hotter, 
the EWs of both ions decrease, although the EW(2600) 
value decreases much faster with increasing temperature 
than EW(2796). 

We now measure Fe II EWs and analyze Fe II line 
profiles to constrain the origin of the Fe II absorption in 
our spectra. 

7.1. Fe II Absorption in TKRS2158 

As noted in N3.ll we observe blueshifted Fe II lines in 
the spectrum of TKRS2158. Velocity offsets are given in 
Table [TJ Velocities are slightly inconsistent among the 
various Fe II transitions; this may be due to poorly sub- 
tracted sky emission, intrinsic Fe II* emission, or other 
sources of noise in the spectrum. In general, Fe II ve- 
locities are slightly lower than the offsets measured for 
Mg II. This may occur if Fe II absorption is weaker at the 
highest velocities, and / or if Mg II emission fills in the 
line profile near systemic velocity and effectively shifts 
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Fig. 11. — Comparison between the outflow absorption strength, 
SFRtjv (left) and M* (right). Each point corresponds to a sin- 
gle coadd. Gray open diamonds plot Wdi f f measurements for the 
TKRS subsamples divided by the 25th and 75th percentile values 
of SFRtjv and M*. Black diamonds show the outflow EW re- 
sults from the decomposition procedure for the same subsamples. 
The red diamond shows the outflow EW (from the decomposition 
procedure) for the "DEEP2-like" subsample selected by M* and 
SFRtjv- Blue circles show results from W09. 2cr upper limits 
are shown when the error on the outflow absorption measurement 
(Wdiff or outflow EW) is twice the central value. 

the line center further to the blue. From Figure [3] we see 
that a few of the Fe II lines extend to the same depth as 
the Mg II lines. We conclude that Fe II traces outflowing 
gas in this object, although we cannot confirm that it 
traces gas in the same velocity range as Mg II. 

7.2. Fe II Absorption in Coadded Spectra 

To analyze Fe II absorption in our galaxy spectra, we 
create a new subsample, selecting objects with z > 0.8 
such that spectral coverage of Fe II AA2586, 2600 is avail- 
able, and further selecting objects with log SFRtjv > 1, 
for a final subsample size of 66 objects. This latter cut 
selects objects with the bluest colors and the strongest 
continua near 2600 A, which minimizes unphysical effects 
due to poorly determined sky levels. Coadding spectra 
of fainter galaxies in the sample yields Fe II absorption 
with unphysical line profiles, and thus are excluded from 
this analysis. 

Figure [TH] shows Fe II and Mg II line profiles for this 
subsample. We note that the Fe II lines are slightly 
deeper than the Mg II lines in the coadd, in con- 
trast to examples fro m QSO absorption line systems 
(jChurchill et aLlfeOOOf ) and stellar atmospheres. This in- 
consistency suggests that a) the depth of Fe II is affected 
by unphysical artifacts resulting from poor sky subtrac- 
tion even in this subsample of the brightest galaxies, b) 
Fe II-absorbing gas has a larger Ct than Mg II, or c) 
the Mg II absorption is being filled in by emission in the 
same transition. 
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Fig. 12. — log7rR^ 2 vs. log SFRtjv for the blue galaxies in 
the TKRS sample. Lines show our subdivisions in logSgpR. All 
galaxies with logEgpri < —1.36 are above the red line; galaxies 
with —1.36 < logXspR < —0.66 are between the red and cyan 
lines; and galaxies with logEgpri > —0.66 are below the cyan line. 
A line of constant logSgpn = — 1 is shown in solid green. The 
dashed green line indicates where galaxies in the sample lie with 
respect to a line of constant logEgp^ = — 1 when the SFRs are 
increased by 0.21 dex. 

We measure the EW in the Fe II lines in the coadded 
spectrum between —300 km s _1 < v < 300 km s _1 , as 
well as on the blue (—300 km s" 1 < v < km s _1 ) 
and red (0 km s _1 < v < 300 km s -1 ) sides of each. 
The total EWs for each line are approximately the same 
(EW(2586) = 1.99 ±0.15 A and EW(2600) = 1.86 ±0.15 
A), indicating that Fe II is completely saturated in this 
coadd. The total EWs of the Mg II absorption lines in the 
same velocity range are significantly larger (EW(2796) = 
2.36 ± 0.10 A), yielding EW(2796)/EW(2600) - 1.2, a 
slightly smaller ratio than expected from stellar atmo- 
spheres and QSO absorption lines. This again may be 
due to emission filling in the Mg II absorption line pro- 
files. We measure a larger EW on the blue side of the 
2600 A line than on the red side (EW(2600, blue) = 
1.12 ± 0.11 A and EW(2600, red) = 0.79 ± 0.10 A). This 
in consistent with a scenario in which the Fe II absorbing 
gas is outflowing. We find that the EW on the blue and 
red sides of the 2586 A line in the coadd are consistent 
within the errors (EW(2586, blue) = 0.99 ± 0.11 A and 
EW(2586, red) = 1.10±0.11 A). We cannot conclude that 
Fe II absorbing gas is outflowing from analysis of this line 
alone; however, the symmetry of the profile does not rule 
out an outflow scenario. The oscillator strength of this 
transition may simply be too lo w for it to obvious l y trace 
blueshifted gas. In contrast, Mart in k. Bouchl (2009) 
find that both the Fe II 2586 and 2600 A lines are satu- 
rated and have similar line profiles in ULIRG outflows. 
See H8. 1.21 for a discussion of the upper limit on N(H) in 
the outflow derived from Fe II EWs. 

We also examine the relative velocity extent of the Fe II 
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Fig. 13. — Sections of the coadded spectra for subsamples divided 
by Esfr around Mg II (black). The coadds have been normalized 
to the level of the continuum surrounding the absorption line s as 
described in the text. Colored lines are the same as in Figure [TOl 
The number of spectra in each subsample is shown in the upper 
right of each panel. 
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Fig. 14. — Outflow absorption strength vs. EgpR. Gray open 
points plot Wdi f f measurements for the subsamples divided by the 
25th and 75th percentile values of EgFR- Black points show the 
outflow EW results from the decomposition procedure for the same 
subsamples. 



and Mg II absorption lines. In Figure [TSJ the Mg II 
2796 absorption extends to larger velocities than the Fe II 
2600 absorption. The Mg II 2796 line has a high velocity 
tail which decreases in strength gradually with increasing 
velocity offset, whereas the Fe II 2600 line has no "tail" 
to high velocities. If the Fe II absorbing gas is indeed 
outflowing, this suggests that the majority of this gas 
does not attain velocities as high as the Mg II absorbing 
gas. 

To quantify this, we measure the relative velocities 
where the profiles of each of these lines reaches a thresh- 
old amount of absorption. We first smooth the coadds, 
and calculate the velocity at which the absorption de- 
creases to 80% and 60% of the continuum level. These 
velocities are marked in Figure [T8l and are closer to the 
center of the line in the case of Fe II, confirming our find- 
ing that Mg II absorption extends to higher velocities. If 
we make the assumption that Fe and Mg have the same 
relative abundances and C/ at all gas velocities, this sug- 
gests that the density of the outflowing clouds decreases 
with increasing veloci ty, as the Fe II column decreases 
with density (see, e.g. JNaravanan et al.H200"8l for ioniza- 
tion modeling of Mg II absorbers). 

8. DISCUSSION 
8.1. Physical Characteristics of Outflows 
8.1.1. Mgll 

Similarly to W09, we use the ratio of the EWs on the 
blue sides of the Mg II lines to estimate a lower limit 
on the column density (iV(Mg II)) in the outflowing gas. 
The oscillator strengths of the lines in the Mg II doublet 
are in the ratio 2:1, and in the optically thin case the 
ratio of the EWs of the lines will also be 2 : 1. As the 
optical depth at line center (to) increases and the lines 
become saturated, the EW ratio approaches 1:1; there- 
fore the EW ratio can be used to constrain To- Once To 
is known, 7V"(M g II) can be calculated using the equation 
(|ST3itzerlll968h 



log N = log 



EW 
log Xf 



log 



2F(r ) 



log C f + 20.053, 



(8) 

adjusted to include the effect of the covering fraction, 
Cf, where A and EW are in A, and N is in atoms cm~ 2 . 
F(t ) is given by: 



F(t ) = / ( 1 _ e ~^ocxp(-a ; )) dx _ 



(9) 



The EW ratio in the Mg II doublet is almost exactly 
equivalent to the "doublet ratio", F (2tq) / F (tq) . After 
this doublet ratio is calculated by taking the ratio of the 
EWs, one may numerically solve for To. This method 
is strictly appropriate only when one absorbing cloud is 
considered. However, it has been found to yield good 
results even when the absorption is caused by a number 
of clouds if the optical depth in the weaker line is To < 5 
(|Jenkinslll98l . 

We measure the EWs for the coadds in the interval 
—384 km s _1 < v < —200 km s _1 for each line in the dou- 
blet; these measurements are listed in Table[5J We choose 
this interval to avoid measuring absorption at the sys- 
temic velocity, which is uncertain to within 200 km s _1 
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Fig. 15. — Normalized SSFR distributions for galaxies in the low-M* (left), middle-M* (middle), and high-M* (right) subsamples. The 
solid histograms show the SSFR, distributions for TKRS galaxies, and the dotted histograms show the SSFR distributions for the W09 
M, -divided subsamples, normalized to have the same maximum values as the corresponding TKRS distributions. 
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Fig. 16.— logSFRuv vs. logM* for the full TKRS sample (in 
black) and a randomly selected portion of the DEEP2 sample (in 
green) from W09. The line shows the cut applied to the TKRS 
sample for the purposes of selecting TKRS galaxies with SFRtjvs 
similar to the SFRtjv values in the W09 sample. 



(see tj2.3p . and to avoid including absorption from the 
red side of the 2796 A line in our 2803 A line measure- 
ment. For the coadd of the entire sample, the EW ratio 
is 1.09 ± 0.16. This yields To, 2803 ~ 3 — oo, where the 
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Fig. 17.— Sections of the high-SFR "DEEP2-like" sample coadd 
around Mg II. The coadd has been normalized to the level of the 
continuum surrounding the absorption lines as des cribed in the 
text. The colored lines are the same as in Figure [Till The number 
of spectra in the sample is shown in the upper right. 



lower limit corresponds to the la upper limit on the EW 
ratio. In the case of the high-M* coadd, the doublet ratio 
is 1.28 ± 0.16, which results in r ,2803 ~ 1-4 — 10. Using 
the equation for N above, and assuming C/ = 1, we find 
log7V(Mg II) > 13.7 for both the full sample and the 
high-M* subsample. Assuming C/ = 0.5 increases log TV 
by 0.3 dex. 
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To estimate a lower limit to iV(H) in the outflow, we 
assume a more conservative value of Cf = 1, and we 
assume AT(Mg) = AT(Mg II); i.e., we do not apply an 
ionization correction. We assume the solar value for the 
abundance of Mg, logMg/H = -4.42, and a factor of -1.2 
dex Mr depletion onto du st, measured in the local ISM 
(|Savage fe; Sembachl 119961 ). This results in the estimate 
logiV(H) > 19.3. We emphasize that this is a very con- 
servative lower limit on the column of outflowing gas be- 
cause of our assumption about Cf, our neglect of ioniza- 
tion corrections, and because our method underestimates 
N for highly saturated Mg II-absorbing velocity compo- 
nents. In addition, the effect of emission from the Mg II 
2796 transition is to reduce the EW measured blueward 
of the Mg II 2803 transition, further reducing the calcu- 
lated log AT (H). These results are nearly an order of mag- 
nitude smaller than the outflow column obtained in W09 
(logiV(H) = 20.1); however, this latter measurement was 
calculated using the EW ratio in the outflow line pro- 
file rather than in the observed profile, and used a wider 
range of velocities (—768 km s _1 < v < km s _1 ). This 
yielded higher EWs in each line and generally higher op- 
tical depths than our la lower limits on r. 

We estimate the mass outflow rate by assuming a spe- 
cific geometry for the outflowing gas. For a thin shell, 
the mass outflow rate is given by 

M « 22 M yr-C /I ^ffl z ,^_^ FT (lO ) 

from W09. The assumption for C/ is unimportant here, 
so long as it matches what was assumed for the calcu- 
lation of N(R) (i.e., the factor of Cf cancels out). We 
have no constraint on the radial extent of the wind from 
our data (R), except that it is likely comparable to the 
size of the galaxies, because C/ is high. For these pur- 
poses, we assume a minimum radius for the shell equal 
to the median half-light radius for the galaxy sample, 
4.1 kpc. We use the velocity at 80% of the continuum in 
the coadd of all galaxies calculated as described in 37.21 
—334 km s _1 . Note that this is not the same velocity 
measurement used in W09, who measured the velocity 
at 50% opacity in the outflow component after decom- 
posing the line profile. The resulting mass outflow rate is 
M > 5 M© yr -1 . As in previous work (e.g.. lMartTnlll999l : 
IRupke et al.ll2005bHMartin|[2005l) . we find an M on the 
same order as the SFR of the sample (~ 1 — 30 M© yr" 1 ). 
We note that M in this sample is lower than in the W09 
sample at z ~ 1.4, and that this is due to lower out- 
flow AT(H) at the lower redshifts. This suggests that the 
galaxies at lower z are less effective in driving outflows 
than at z ~ 1.4. However, we reiterate that the outflow 
AT(H) in both this work and in W09 may be substantially 
underestimated. 

8.1.2. Fell 

Because the profiles of the two Fe II lines do not 
have an asymmetric blue wing, we hypothesize that ei- 
ther these lines are too weak to trace outflows at ve- 
locities > 300 km s _1 , or alternatively that Fe II sim- 
ply is not present in outflowing gas at these velocities. 
However, assuming the absorption on the blue sides of 
these lines at lower relative velocities is due to outflow- 
ing gas, we again use an EW ratio, this time in the ve- 



locity range —300 km s" 1 < v < —200 km s" 1 , to cal- 
culate a lower limit on the outflow column in the Fe II 
transition. The EWs of the two lines in this velocity 
range for the coadd we use to analyze Fe II are both 
EW = 0.14 ± 0.06, yielding an observed EW ratio of 
1.0 ± 0.7. The two Fe II lines have oscillator strengths 
in the ratio 3.475:1. We numerically solve the equa- 
tion F(3.475t )/F(t ) = 1.0 ± 0.7 for r , 2 586 and ob- 
tain the limit To,2586 }t 1-5. Using logA2586/ = 2.252 
(|Mortonll2003l ) and EW b k,c(2586) = 0.08-0.2 A, we find 
logA r (Fe II) > 13.5 — 13.9. Assuming solar abundance, 
no ionization correction and a factor of 1/10 depletion 
onto dust, we obtain logAT(H) ~ 19.0 — 19.4, in good 
agreement with our result from Mg II. In this case, anal- 
ysis of the Fe II profiles does not yield a more stringent 
constraint on the outflow column than analysis of Mg II. 

We may compare column densities measured in differ- 
ent relative velocity ranges for these lines if we assume 
that the distribution of N in individual clouds changes 
smoothly with changing velocity and if we assume that 
the Cf of the outflow is constant at all velocities. We 
calculate a 3er upper limit of ~ 0.3 A on the EW of ab- 
sorption at velocities —500 km s" 1 < v < —300 km s _1 
with respect to both Fe II absorption lines. We use this 
measurement to set an upper limit on the outflow col- 
umn at these velocities by assuming that this gas must 
have an optical depth lower than the limit determined 
above for the range —300 km s _1 < v < —200 km s _1 . 
Again we use the equation for log N given in the pre- 
vious section and apply it to the case of the 2600 A 
line, with r n . 2 finn = 3.475t , 2 586 = 5.2 and logA 26 oo/ = 
2.793 (Morton 2003). We find logA^(Fe II) < 13.8 and 
logAT(H) < 19.3, using the same assumptions for abun- 
dance, ionization correction and dust depletion as above. 
Note that this result requires that Ct — 1 at all velocities; 
this assump tion does not hold for th e ULIRG outflows 
discussed in lMartin fe Bouchel (j2009f l. 

8.2. Outflow Absorption Dependence on SFR and M* 

We have shown that outflow absorption strength in- 
creases with SFRuv and M„ in the TKRS sample. This 
statement addresses the behavior of outflow absorption 
in the middle- and high-SFRuv and M» subsamples 
only (with S/N > 6 pix -1 ); in the subsamples with 
the lowest values of these parameters, the coadds have 
S/N ~ 3 pix -1 and only poorly constrain outflow prop- 
erties. It does appear that the lowest-M* subsample has 
a large enough outflow EW to be significant (> 1 A; 
see discussion in §4.3|k however, the low measurement of 
Wdiff is consistent with no outflow. 

Interestingly, both the middle-SFRuv an d M* sub- 
samples have Wdiff and outflow EW measurements < 
0.4 A, consistent with no outflow. These coadds have 
S/N ~ 6 pix -1 ; in the simulated spectra discussed in 
34.31 saturated ("strong") outflows with v > 100 km s _1 
and br> > 150 km s _1 yield values of outflow EW greater 
than 0.85 A with 95% confidence. Certainly weaker 
outflows with smaller bo, velocities < 100 km s -1 , and 
with column densities such that the total EW(Mg II 
2796) < 3.5 A, would not be detected in these coadds, 
and such outflows could exist in these galaxies. How- 
ever, we find that the middle subsamples of galaxies with 
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0.555 < logSFRuv < 0.942 and 9.86 < logM* < 10.49 
at z ~ 1 are not driving strong outflows. 

Galaxies with M* and SFRyv above these ranges show 
clear evidence for outflow. This includes galaxies in our 
subsample specifically selected to have similar SFRuv 
and M* to the W09 sample at z ~ 1.4 (the "DEEP2- 
like" subsample). It also includes the galaxies in the W09 
sample, which all have logSFRuv > 1- Our middle-M* 
subsample overlaps with the W09 sample in M* but falls 
below it in SFRuv, suggesting that outflow absorption 
strength is more closely associated with SFR than with 
M* at z = 1 — 1.4. We also note that outflow absorp- 
tion strength in general does not increase with increas- 
ing specific SFR (SSFR). Figure QJ shows the SSFR-M* 
distribution of the TKRS and W09 M* -divided samples 
with black and blue diamonds, respectively. Although 
the middle-M, TKRS subsample has higher SSFR than 
the high-M* TKRS subsample, it exhibits no outflow sig- 
nature (see Figures [TU] and [TTJ) . 

It is interesting to consider whether outflowing gas 
will escape from these high-M* and high-SFRuv galax- 
ies. The escape velocity of galaxies can be estimated 
from the width of the [O II] emission line as in W09. 
The [O II] linewidths (a) of TKRS galaxies were mea- 
sured in iWeiner et alj (|2006f h We find that the median 
a for both our high-M* and high-SFRuv subsamples is 
cr([0 II]) = 80 km s _1 , and use the relation from W09: 



escape 



5-6cr([0 II]). 



(11) 



This yields an escape velocity ~ 480 km s" 1 . From the 
outflow components plotted in green in the bottom pan- 
els of Figure [TUJ we see that in the high-SFRuv sub- 
sample, which is most like the W09 sample in terms 
of SFR, there appears to be a tail of Mg II absorption 
beyond the escape velocity. In the other subsamples, 
there is no compelling evidence of gas beyond the es- 
cape velocity. We conclude that is it unlikely that the 
gas detected in absorption is able to escape from most 
of the galaxies in our sample. This differs from the re- 
sults of W09, who find that the absorption line profile 
extends well beyond the median escape velocity for each 
of their M*-divided subsamples. The profile extends to 
particularly high velocities in the highest-M* galaxies; 
the velocity where the profile is at 90% of the contin- 
uum level exceeds the escape velocities for > 84% of the 
subsample with logM* > 10.45. We cannot conclude 
from this analysis, however, that cool Mg II- absorbing 
gas does n ot escape from the bulk o f the TKRS sample 
at z ~ 1. iMartin fe Bouchel (|2009f) found that in low- 
ionization ULIRG outflows, the absorption is saturated 
at all velocities, while the C/ of the gas decreases with 
increasing outflow velocity. This implies that the veloc- 
ity extent of the detected absorption profile is controlled 
by the decreasing gas C/ and the S/N of the data, rather 
than decreasing N at high velocities, and that there may 
be very low Cf clouds moving at well above the escape 
velocity. Higher S/N coadds would be needed to con- 
strain the absorption depth at the escape velocity more 
precisely. Hotter, and likely more diffuse, gas kinematics 
are not constrained by our observations or those of W09. 

8.3. Trends in Outflow Absorption with Ssfr 

As discussed in [Hcckman (2002f) noted that lo- 
cal galaxies which drive winds strong enough to be de- 
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Fig. 18. — Sections of the coadd of galaxy spectra used in the 
Fe II analysis with logSFRuv > 1. The black spectrum shows 
the region surrounding Fe II 2586, 2600, and the blue spectrum 
shows the region surrounding Mg II 2796, 2803. Coadds have been 
normalized to the level of the continuum surrounding the absorp- 
tion transitions, indicated with the horizontal lines. The section 
surrounding Mg II has been offset by +1 in flux. The velocity for 
the black spectrum is measured relative to the rest velocity of the 
2600 A line; for the blue spectrum it is measured relative to the 
rest velocity of the 2796 A line. Vertical red dotted lines show rest 
velocities of all lines of interest. Vertical blue and black marks in- 
dicate the velocities at which the absorption decreases to 80% and 
60% of the continuum level for the Mg II 2796 and Fe II 2600 lines, 
respectively. 

tectable in absorption universally satisfy the criterion 
Esfr > 0.1 M Q yr" 1 kpc~ 2 . The galaxies in the W09 
sample and LBGs at z ~ 2 — 3 (Sha plev et alj 120031 : 
iPettini eT al. 2002) also meet this criterion. Further evi- 
dence for a threshold Xsfr is provided by studies of out- 
flow remnants in emission. For instance, in a study of 
the relationship between the sizes of radio halos created 
by energy input from star formation and t he sizes of star- 
formin g regions in local spiral galaxies, iDahlem et al.l 
(2006) found that radio halo size is linearly correlated 
with the size of star-forming regions. They additionally 
found that there is a threshold energy input rate per 
unit surface area above which radio halos form, and that 
this threshold depends on the mass surface density of the 
galaxy. 

There are two different physical mechanisms that 
may be responsible for driving these winds, which re- 
sult in different dependencies of outflow properties on 
SgFR- One mechanism is the thermalization of en- 
ergy from supernova ejecta, which generates a hot wind 
that e ntrains cold clouds i n an o utflow via ram pressur e 
(e.g.. IChevalier fe Cleed [19851: IHeckman et ail IWXb . 
iMcKee fc Ostrikerl(|1977t ) suggested that a threshold spa- 
tial density of supernovae must be achieved in order 
for this mechanism to operate. Regions with super- 
novae densities above such a threshold create and main- 
tain a superbubble filled with hot gas with a cooling 
timescale much longe r than the bubble expansion time 
(jHeckman et al.fl 9 901) . Further, at higher supernova rate 
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Fig. 19. — Diamonds show SSFR vs. M* for the W09 and TKRS subsamplcs divided by M» (blue and black, respectively). The point 
size is scaled to the square root of the outflow EW for each subsample. Error bars show the ±34th-percentile values of SSFR and M« for 
each subsample. The open diamond indicates a subsample for which analysis of the low-S/N coadd provides only an upper limit on the 
outflow EW. Dotted lines show fits to the SSFR-M, main sequence derived in Noeskc et al. (2007a) at z t, s = 0.35, 1 and 1.4 (bottom to 
top). Solid lines show the paths of r model galaxies through the diagram with time for logM;, = 10.5, 10.75, 11.0 and 11.25 (left to right). 
The open circles along these lines indicate the location of individual models with different formation redshifts (zf) at z (, s = 0.35, 1.0 and 
1.4. Each color corresponds to a model with a different formation redshift, and thus generally appears 3 times along each solid line (i.e., 
there is one circle of each color for each observed redshift). See text for more details. 



per unit volume (or at lower ambient gas density) , radia- 
tive losses in the bubble are reduced, and the efficiency of 
thermalization of the e.jecta is enhan ced (Tflcckma n et al.l 
119901: iStrickland k. Heckmanl 12009ft . This spatial su- 
pernova rate density is closely tied to the ob servable 
quant ity, E sfr.- Additionally, a s described by iMartinl 
(2005) and Murray et all ((2005), in this energy-driven 
paradigm, the terminal velocity of cold clouds in the out- 
flow is proportional to Rq ' 5 , where Rq is the size of the 
starburst. 

The other mechanism that may be important in driv- 
ing winds is momentum from radiation pressure or cos- 
mic ray pressure generated by a starb urst or AGN 
(|Murrav et al.l l2005t iSocrates et~ai1 [2008) . In the case 
of radiation-driven outflow, the dependence of outflow 
velocity on starburst radius for an optically thick wind 



is slightly weaker, (In l/i?o) ' 5 , while the criterion that 
must be met to drive a wind depends on the masses of 
the luminous star clusters driving the outflow. The star 
cluster mass, like spatial supernova density, is closely re- 
lated to Esfr (see N. Murray & B. Menard, in prepara- 
tion, for further discussion). Thus, while both proposed 
mechanisms require that a threshold Esfr be satisfied 
to take effect, the empirical relationship between out- 
flow velocity and starburst radius may eventually help 
constrain the relative contributions of these two mecha- 
nisms for driving outflows, although these relationships 
are complicated by additional factors (e.g., viewing an- 
gle, outflow morphology and Ct, detection technique and 
sensitivity). 

Our sample suggests that outflow absorption strength 
increases with Esfr- As in £18.21 although we measure 
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relatively high values of Wdiff and outflow EW for our 
lowest-EgpR coadd (with logSsFR < —1-36), because 
this spectrum has low S/N, we consider the detection 
of outflow tentative, and it must be confirmed in future 
analyses. In the middle and high-SsFR subsamples, the 
difference between outflow absorption measurements is 
significant only in the values of Wdiff, the outflow EW 
values are consistent within the errors. It is interesting, 
though at the limit of our data, to test whether the local 
criterion on Ssfr for the driving of outflows also ap- 
plies at z ~ 1. By accounting for a possible systematic 
underestimate of SFR by 0.21 dex (see §5.4|) . we have 
constructed a subsample with a Ssfr that is likely to be 
under the local threshold in the mean. This coadd has 
S/N ~ 4.7 pix" 1 , and yields an outflow EW ~ 1 A. This 
value is su gges tive of outflow, given the results of our 
analysis in §4.3[ although the outflow absorption strength 
is not as large as in the high-M* and high-SFRuv sub- 
samples. However, given the minimum estimated disper- 
sion in our SFRs of 0.43 dex, this subsample may include 
some galaxies with Ssfr ^ 0.1 M yr _1 kpc~ 2 . 

Additionally, we have measured only a flux-averaged 
SsfRj with galaxy sizes parametrized by the half-light 
radius in the UV. Starburst radii in the local Uni- 
verse are often constrained using half-light radii mea- 
sured from UV or Ha imaging; however, local and dis- 
tant starbursts have been found to have distinct mor- 
phologies which may make measurements of half-light 
radii less appropriate for this purpose at z *~ 1. For 
instance, local LIRGs commonly have centrally concen- 
trated star formation (Sanders & Mirabel! Il996f ). while 
LIRGs at z ~ 1 have more widely distributed star- 
forming knots (|Melbourne et alj 120081 ). The half-light 
radius may therefore not be an accurate characteriza- 
tion of starburst radius in the distant Universe, where 
it reflects the sizes of the distributions of star-forming 
knots in the objects, rather than the total surface area 
of the knots themselves. Although we have not detected 
a strong dependence of outflow absorption strength on 
Ssfr, we have by no means ruled it out, and future stud- 
ies with higher S/N spectra should include more detailed 
analyses of Ssfr in distant galaxies. 

8.4. Redshift Evolution of Outflow 
8.4.1. Following the Sample Evolution 

We have identified evolutionary trends in outflow ab- 
sorption strength which are dependent on galaxy M*. 
Specifically, we have found that at logM* > 10.49, out- 
flows are of comparable strength at z = 1 and 1.4. At 
lower M*, however, there is a significant decrease in out- 
flow absorption strength with decreasing redshift along 
with a decrease in the SFR in the same galaxies. Stated 
in another way, galaxies with similar SFRs at z = 1 
and 1.4 have similarly large outflow absorption strengths, 
while galaxies with logSFR uv < 0.942 (< 8.7 M Q yr" 1 ) 
do not appear to drive strong outflows. While these 
conclusions apply to galaxies in specific and SFRuv 
ranges at different redshifts, they do not speak to the evo- 
lution of the outflow properties in a given galaxy through 
time. We now attempt to associate the TKRS galaxies 
with their progenitors in the W09 sample at z ~ 1.4 by 
creating so-called "r models" to parametrize the evolu- 
tion of SFR and M* in a given galaxy. We are largely 



motivated by the success of r models in reproducing the 
distribution of galax ies in SSFR-M* space out to z ~ 1 
(jNoeske et al.ll2007al) . These models are too simplistic to 
fully describe the evolution of an individual galaxy, which 
may experience starbursts, mergers, and other distur- 
bances that define its star formation history. However, r 
models were intended to describe the broad evolution of 
populations of galaxies, and are appropriate here, as we 
are interested in the evolution of large galaxy subsam- 
ples. 

In the context of this model, the SFR of a galaxy of 
total baryonic mass Mb at redshift z is given by 



SFR(M b , z) = SFR(z/) exp 



(12) 



where Zf is the formation redshift of the galaxy, T — 
t(z) —t(zf) is the gal axy age, and r is the e-folding time. 
INoeske et "all (|2007aD found that models with 



r(M b 



c a M£ 



and 



I + z f (M b ) = cpM; 



(13) 



(14) 



with the constants 10 20 4 < c a < 10 20 7 , 10~ 2 - 5 < eg < 
-1, and (3 = 0.3 successfully reproduce the 



10- 1 - 7 , a = 

SSFR-M* "main sequence" as well as the spread around 
this sequence for galaxies from the AEGIS survey at 
0.2 < z < 1.1. The changing dependence of Zf on M b is 
necessary to account for the high SSFRs of galaxies at 
low M* . Although these models have not been tested out 
to z ~ 1.4, where the AEGIS sample is incomplete even 
at the highest M* , we assume that they describe the full 
locus of star-forming galaxies at this redshift. 

Figure [TO] shows the location in SSFR-M* parame- 
ter space of the W09 M*-divided subsamples (blue dia- 
monds) and the TKRS M^-divided subsamples (black di- 
amonds) . The point size increases with the square root of 
the outflow EW for each subsample. The open diamond 
indicates the location of the lowest-M* TKRS subsample, 
which has insufficient S/N to yield an outflow detection 
(at 1.2 A) with high confidence. The error bars indicate 
the spread in galaxy properties in each subsample at the 
±34th percentiles. The loci of the SSFR-M* "main se- 
quence" at observed redshifts z b s — 0.35, 1.0 and 1.4 are 
shown with nearly horizontal dotted lines from bottom to 
top. These loci d escribe the locations of galaxies in the 
AEGIS sample of INoeske et"aTl (|2007al ) at z < 1.1. The 
four solid lines show the path of r model galaxies of var- 
ious M b (with models for logM 6 = 10.5, 10.75, 11.0 and 
11.25 plotted from left to right). The circles along the 
single r model paths mark the locations of model galaxies 
with different formation redshifts (z/) at observed red- 
shifts z bs = 0.35, 1.0 and 1.4. The Figure legend gives 
the appropriate Zf for each color; for instance, purple cir- 
cles are for the latest Zf = 1.8 models, and green circles 
are for the earliest zj = 5.0 models. 

The W09/DEEP2 galaxies lie slightly above the SSFR- 
M* main sequence at z ~ 1.4. This is because the DEEP2 
i?-band selection limits the sample to the highest-SFR 
galaxies at a given M*; the effect is most severe at 
z = 1.4, the redshift limit of the survey. (As noted previ- 
ously, the model for the main sequence at z ~ 1.4 has not 
been observationally tested and is an extrapolation from 
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lower redshift models.) The TKRS galaxies, on the other 
hand, lie below the main sequence at z ~ 1. This discrep- 
ancy m ay be due to a number of factors. iNoeske et alj 
(2007b) used a combination of Lir and emission line 
luminosities to constrain the total extinction-corrected 
SFRs in their sample. We have shown that our UV- 
based SFRs are systematically lower (by 0.21 dex) than 
IR-based SFRs for galaxies with IR detections; thus there 
may be a s ystematic offset betw een the SFRs in this 
work and in INoeske et alJ (|2007bl ). In addition, we may 
have failed to remove some red, non-star-forming galax- 
ies from the highest-M* subsample, which would reduce 
the median SFR in this bin. While the offset between the 
main sequence and our sample is as large as ~ 0.25 dex 
in SSFR for the highest-M» subsample, this discrepancy 
does not significantly affect our conclusions, discussed 
below. 

8.4.2. Outflows in the Sample Progenitors 

We first consider the highest-M* TKRS subsample. 
The median galaxy in this sample has Zf ~ 5; i.e., it 
lies near the green model along the log Mb = 11.0 evo- 
lutionary track at z a bs = 1- I n the case that we are 
underestimating the sample SFR by ~ 0.2 dex, the me- 
dian galaxy would lie near the cyan model with Zf ~ 3. 
The highest-M* W09 subsample appears to include the 
progenitor of this median galaxy (the green or cyan mod- 
els at z {,s = 1.4), as well as the highest SSFR galaxies in 
the subsample (with later zf). It appears that outflows 
persist at similar strengths in high M* galaxies between 
z = 1.4 and z = 1. Even massive galaxies with early 
Zf > 5 and very low SSFRs are apparently still driving 
outflows at z ~ 1. 

We next turn to the lower-M* subsamples. The high- 
SSFR galaxies in the middle-M, TKRS subsample lie 
close to the logAf;, = 10.75 evolutionary track. There 
are W09 galaxies that lie along this track (in the lowest- 
or middle-M* W09 subsample), which have SSFRs that 
may be consistent with progenitors of these high-SSFR 
TKRS galaxies, particularly if we consider that our SSFR 
values may be systematically low. That is, a r model 
galaxy with log Mb = 10.75 and Zf = 2 (shown in blue) 
may be included in the middle-M* W09 subsample when 
it is observed at z b s = 1.4 and passes through the loca- 
tion of the middle-M* TKRS subsample at z j, a = 1. In 
the case that we are underestimating the SSFR of this 
subsample, the middle-M» galaxies at z ~ 1 will lie close 
to the purple model with Zf = 1.8; the progenitors of 
this model pass through the lowest-M* W09 subsample 
at z b s = 1.4. The galaxies with the lower SSFR val- 
ues in this middle-M* TKRS subsample, however, have 
no counterpart in W09. We cannot constrain from this 
analysis what fraction of the galaxies in this subsample 
have progenitors in the W09 sample. However, the W09 
galaxies were found to ubiquitously host outflows, while 
the middle-M* TKRS subsample does not host strong 
outflows. It is unclear whether the outflows at z ~ 1.4 
have been shut off by z ~ 1, or if the lack of outflows in 
lower-SSFR galaxies is diluting the outflow signature. 

The galaxies in the lowest-M» TKRS subsample may 
have a few progenitors in the low-M* end of the lowest- 
M* W09 subsample; i.e., there are a few W09 galaxies 
that lie between the log Mb = 10.5 and 10.75 evolution- 
ary tracks. However, we do not comment on the presence 



or absence of outflows in these galaxies at z ~ 1, as we 
lack robust outflow measurements for this subsample due 
to the low S/N ~ 2.3 pix -1 of the coadd. 

In summary, we find that outflows persist in galaxies 
with logM, > 10.5 between z ~ 1.4 and 1, while in 
our middle-M* subsample, outflows may have shut down 
during this time. We lack the S/N to comment on the 
evolution of outflows in the lowest-M* TKRS galaxies. 

8.4.3. The Sample Descendants 

ISato et al . (2009) investigated outflows traced by Na I 
absorption in AEGIS galaxies at 0.11 < z < 0.54. Keep- 
ing in mind the caveat that Na I likely traces cooler, 
higher density gas than Mg II, we may search f or the 
descendants of TKRS galaxies in the ISato et all (I2009T) 
sampl e and compare outflow properties. ISato et all 
(2009) are able to measure outflow properties only in 
galaxies with the highest S/N near the Na I doublet; 
these are th e highest- S/N spectra in the sample o f 
INoeske et all (|2007al lH). As a result, ISato et all (|2009f) 
are complete only down to log M t > 10.75, and only 
at the lowest redshifts (z < 0.3). Among these galax- 
ies, the vast majority have Na I absorption at the galaxy 
systemic velocity, or redshifted with respect to systemic 
velocity. At z > 0.3, the galaxies with sufficient S/N to 
measure Na I kinematics have a mix of outflows, inflows, 
and absorption at systemic velocity. These galaxies are 
massive enough to be descendants of the galaxies in the 
high-M* W09 and TKRS subsamples, and yet do not 
appear to be driving Na I outflows ubiquitously, particu- 
larly at the lowest redshifts. It may be that either these 
galaxies are driving outflows of slightly warmer gas with 
lower density than can be probed by Na I, or that their 
outflows have been halted, possibly because the ISM has 
been removed. These galaxies have lower SFRs than 
massive galaxies at z ~ 1, consistent with the idea that 
their cool gas supply has been reduced or exhausted. 

In galaxies with logM, < 10.75, many of which may 
be descendants of the lower-M* TKRS subsamples, there 
is again a range in Na I kinematics. Most of the galax- 
ies with kinematic information available are not driving 
outflows; however, many more have gas kinematics which 
have not been measured. Thus the outflow properties of 
the descendants of the lower-M» TKRS galaxies remain 
unconstrained. 

Many of the galaxies in the TKRS sample, like those in 
the W09 sample, have log M* < 10.7; i.e., they are be low 
the "quenching mass" derived bv lBundv et al.l (|2006f ) for 
0.4 < z < 0.7. This quenching mass describes the mass 
at which star formation is suppressed, and decreases from 
z = 1.4 to 0.4. These lower-M* galaxies are not driving 
strong winds for the most part, and because they are 
well below the quenching limit, are likely to remain blue, 
star-forming galaxies today. The TKRS galaxies with 
M» above this limit will likely evolve to the red sequence 
by z = 0. 

9. CONCLUSIONS 

We have an alyzed the spectra of 468 galaxies in the 
TKRS survey IjWirth et al]l2004fl of the GOODS-N field 
at 0.7 < z < 1.5 and identified cool outflowing gas 
traced by Mg II absorption in both individual galaxy 
spectra and in coadded spectra. This is the first report 
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of cool outflows in this redshift range. We find that the 
most massive, highest-SFR galaxies in the sample have 
a strong outflow absorption signature, while less massive 
(logM* < 10.5, logSFRuv < 0.94) star-forming galaxies 
have no outflows or outflows with v < 100 km s^ 1 and 
with bo < 150 km s — 1 . These same low-M* galaxies have 
higher specific SFRs (SSFR) than the high-M» subsam- 
ple, but lower absolute SFRs, suggesting that outflow 
absorption strength is more closely correlated with the 
latter. 

We find that the outflowing gas has a column density 
log iV(H) > 19.3, resulting in a mass outflow rate of order 
the SFR of our sample. This column is slightly smaller 
than the outflow column densities observed at z ~ 1.4 in 
the DEEP 2 survey (W09) and in local IR luminous galax- 
ies (e.g., lRupke"et^[2005af) . Our comparison of the 
Fe II and Mg II line profiles in coadded spectra of high- 
SFR galaxies suggests that Mg II emission is a common 
feature of galaxies in this sample. Although it remains 
unclear whether the Fe II absorption in this sample traces 
outflowing gas, by making the assumption that Fe II is 
indeed present in the same outflowing clouds that are 
traced by Mg II and that the C/ is independent of out- 
flow velocity, we use the lack of Fe II absorption blueward 
of —300 km s _1 to place an upper limit on the outflow 
column density at these velocities of log iV(H) ~ 19.3. 

We also investigate, for the first time, trends in outflow 
absorption strength with SFR surface density (Ssfr) a t 
z ~ 1. We find tentative evidence for a weak trend (at 
~ la significance) of increasing outflow absorption with 
increasing Esfr at logSsFR > —1-36. We cannot rule 
out outflows in galaxies with EgpR, below the local thresh- 
old for driving winds, Esfr = 0.1 Mq yr" 1 kpc~ 2 , and 
in fact see some evidence for outflows in these objects. 
Future work in this area should focus on developing a 
set of measurements of the surface area over which star 
formation occurs, carefully calibrated total SFRs, and 
outflow velocities in individual galaxies. 

Finally, we examine the temporal evolution of outflows. 
We find that galaxies with similar SFRs at z ~ 1.4 and 
z ~ 1 have similarly strong outflow absorption. To exam- 
ine the evolution of outflows in specific galaxies through 
time, we invoke the r model star formation histories of 
INoeske et"all (|2007al) . who note that the decline in SFR 
at a given M* with time is consistent with being driven by 
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gas exhaustion. We find that outflows persist in massive 
(logM»/M Q > 10.5) galaxies between z ~ 1.4 and z ~ 1, 
in spite of a decreasing SFR between these t wo epochs. 
Meas urements of Na I kinematics at z ~ 0.2 (Sat o et al.l 
2009) suggest that these massive galaxies, whose SFRs 
will continue to decline over time, will stop driving winds 
at low redshift. We speculate that this is because the 
cool gas supply in the ISM will be exhausted by z ~ 0.2. 
Many (but not all) of these galaxies will become red and 
dead ellipticals by this time. We measure weak or no 
outflows (with outflow v < 100 km s _1 ) in galaxies with 
lower masses at z ~ 1. Based on our comparison of the 
SSFR and M* of our sample to those of r model galaxies 
with exponentially declining star formation histories, we 
find that the progenitors of some of these lower-M* galax- 
ies are likely contained in the DEEP2 z ~ 1.4 compar- 
ison sample and are driving strong outflows at z ~ 1.4. 
These objects are likely to remain blue cloud galaxies 
until z ~ 0. It appears that the outflow strength or 
the frequency of outflows in this population has declined 
from z ~ 1.4 to z ~ 1, as has the SFR. It is undeter- 
mined whether they will ever drive strong outflows as 
they evolve from z ~ 1 to today. 
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TABLE 1 
Measurements of UV 

ABSORPTION LINE VELOCITY 
OFFSETS IN TKR.S2158 



Transition Ad 

km s~ 1 



Fc n 2344 
Fc ii 2374 
Fo ii 2382 
Fc ii 2586 
Fc ii 2600 
Mg ii 2796 
Mg ii 2803 



-262 ± 29 
-281 ± 42 
-415 ± 36 
-287 ± 40 
-229 ± 62 
-314 ± 19 
-257 ± 30 



Note. — Systemic 
velocities are determined 
from the velocity of the 
[O ii] emission line dou- 
blet. 
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EW Measurements of Coadded Spectra 
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Note. — Upper limits are given at the 2a level. No doublet ratio is reported if either of the EW 
measurements is an upper limit. 

a Measured in velocity interval —384 km s _1 < v < —200 km s — 1 relative to the systemic velocity of each 
transition. 



